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fir  AfSTRACT  (MonmomiOOwonM) 

The  characterization  of  reverberation,  from  the  surface,  the  bottom,  and/or  the 
volume,  is  important  in  advancing  the  understanding  of  the  mechanisms  involved  in 
underwater  sound  scattering.  Numerous  experiments  and  much  theoretical  work  have 
been  done  in  the  area  of  monostatic  reverberation,  but  the  area  of  bistatic 
reverberation  has  been  much  less  researched.  This  certainly  is  a  result  of  the  much 
more  complex  geometries  involved  in  bistatic  reverberation. 

This  dissertation  studied  the  statistics  of  bistatic  surface  reverberation  data 
from  the  FLIP  experiment  conducted  in  January  of  1992.  The  data  were  verified  as 
homogeneous  and  normally  distributed  using  the  Kolmogorov-Smimov  two-sample 
test  with  a  confidence  interval,  a,  equal  to  0.1.  using  a  technique  by  which 
ensembles  from  different  times  were  combined,  meaningful  deviations  from  a  normal 
distribution  verr  observed  at  the  highest  wind  speed  of  7.2  m/s. 
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Abstract 


The  characterization  of  reverberation,  from  the  surface,  the  bottom,  and/or 
the  volume,  is  important  in  advancing  the  understanding  of  the  mechanisms  in¬ 
volved  in  underwater  sound  scattering.  Numerous  experiments  and  much  theoret¬ 
ical  work  have  been  done  in  the  area  of  monostatic  reverberation,  but  the  area  of 
bistatic  reverberation  has  been  much  less  researched.  This  certainly  is  a  result  of 
the  much  more  complex  geometries  involved  in  bistatic  reverberation. 

This  dissertation  studied  the  statistics  of  bistatic  surface  reverberation  data 
from  the  FLIP  experiment  conducted  in  January  of  1992.  The  data  were  verified  as 
homogeneous  and  normally  distributed  using  the  Kolmogorov-Smirnov  two-sample 
test  with  a  confidence  interval,  or,  equal  to  0.1.  Using  a  technique  by  which  en¬ 
sembles  from  different  times  were  combined,  meaningful  deviations  from  a  normal 
distribution  were  observed  at  the  highest  wind  speed  of  7.2  m/s. 

The  bistatic  surface  scattering  strengths  were  calculated  from  the  data  and 
compared  with  the  predictions  of  a  theory  developed  by  S.T.  McDaniel  for  the 
prediction  of  bistatic  surface  scattering  strengths.  The  comparison  showed  good 
agreement  with  a  rms  deviation  of  about  3  dB. 
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Chapter  1 


Introduction 


1.1  Motivation 

The  characterization  of  underwater  reverberation,  from  the  surface,  the 
bottom,  and/or  the  volume,  is  important  in  advancing  the  understanding  of  the 
mechanisms  involved  in  underwater  sound  scattering.  Numerous  experiments  have 
been  conducted  to  understand  the  effects  of  the  environment,  e.g.,  wind  speed, 
wave  height,  grazing  angle  and  bottom  roughness,  on  reverberation.  Theoretical 
predictions  or  models  for  reverberation  have  been  developed,  some  based  entirely 
upon  theoretical  work  and  some  using  empirical  methods.  This  thesis  establishes 
the  characteristics  of  bistatic  ocean  surface  reverberation  as  a  function  of  wind 
speed  and  compares  the  mean  bistatic  surface  scattering  strengths  to  one  predictive 
model. 

Reverberation  may  be  classified  according  to  its  source:  surface,  bottom, 
volume,  or  some  combination.  In  addition,  reverberation  may  be  classified  as  ei¬ 
ther  monostatic  (source  and  receiver  collocated)  or  bistatic  (source  and  receiver 
at  different  locations).  To  date,  boundary  reverberation  research  has  focused  on 
monostatic  geometries  with  few  exceptions  [1,  2,  3].  The  data  analysis  in  this  thesis 
focuses  on  the  statistical  characteristics  and  mean  scattering  strengths  of  bistatic 


surface  scatter  as  a  function  of  geometry  and  wind  speed. 

From  a  geometric  standpoint,  monostatic  boundary  scattering  strength  is 
a  function  of  one  parameter,  the  grazing  angle,  denoted  by  v’mc-  It  may  also  be  a 
function  of  wind  direction  2ls  suggested  by  Zornig  [1].  Bistatic  boundary  reverbera¬ 
tion  also  depends  upon  the  bistatic  angle,  denoted  by  <t>,  and  the  scattered  grazing 
angle,  denoted  by  rl>,cat-  There  are,  of  course,  additional  parameters  necessary  to 
describe  dependence  upon  the  environment.  The  additional  geometrical  param¬ 
eters  present  some  challenges  in  displaying  mean  bistatic  scattering  strength,  as 
discussed  in  Section  5.4. 

Characterizing  reverberation  usually  means  describing  the  mean  scattering 
strength,  5,  as  a  function  of  the  geometrical  and  environmental  parameters.  A 
more  complete  characterization  is  contained  in  the  cumulative  distribution  func¬ 
tion,  cd/,  or  probability  density  function,  pdf,  of  S.  The  mean  and  the  variance 
can  be  calculated  from  the  distribution,  and  more  importantly,  comparison  with 
theoretical  distributions  can  be  made.  The  statistical  analysis  of  the  bistatic  sur¬ 
face  scatter  data  implemented  in  this  thesis  used  the  cdfoi  the  probability  of  false 
alarm,  which  is  defined  as  one  minus  the  cdf  'xn  the  analysis  of  reverberation 
data. 

This  thesis  presents  an  analysis  of  the  statistical  characteristics  of  bistatic 
surface  scattering  strength  from  ocean  surface  scatter  data  taken  in  January  of  1992. 
The  experiment  was  conducted  by  the  Applied  Physics  Laboratory  at  the  University 
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of  Washington  (APL-UW)  using  the  RV  FLIP,  Floating  Instrument  Platform. 
The  main  purpose  of  the  experiment  was  to  investigate  the  relationship  between 
the  ocean  sub-surface  bubble  layer  and  surface  forward  scatter  and  backscatter. 
Forward  and  backscatter  data  were  taken  between  10  and  80  KHz.  Out  of  plane 
scatter  was  also  investigated  over  the  range  of  frequencies,  of  which  data  at  30  and 
40  KHz  could  be  analyzed  for  bistatic  scattering  strength. 

1.2  Original  Contributions 

The  qu2dity  and  uniqueness  of  the  data  produced  a  number  of  new  results. 
The  original  contributions  presented  in  this  thesis  are  as  follows.  These  contribu¬ 
tions  are  discussed  in  more  detail  in  Chapters  4,  5  and  6. 

•  The  ocean  surface  scatter  data  were  histogrammed  and  their  cumulative  distri¬ 
bution  function  compared  with  a  Rayleigh  distribution  using  the  probability  of 
false  alarm,  Pfa-  Theory  predicts  that  the  dominant  surface  scattering  mech¬ 
anism  is  the  slope  of  surface  for  higher  grazing  angles  and  isotropic  resonant 
bubble  scatter  at  lower  grazing  angles  and  higher  wind  speeds.  The  magni¬ 
tude  of  the  complex  envelope  of  the  ocean  surface  scatter  data  was  found  to 
be  Rayleigh  for  both  mechanisms  and  for  all  but  the  highest  wind  speed  (7.2 
m/s).  This  result  shows  that  the  use  of  the  central  limit  theorem  to  describe 
the  statistical  characteristics  of  bistatic  surface  scatter  is  valid  over  the  geome¬ 
tries  and  wind  speeds,  except  7.2  m/s,  studied  here.  It  also  shows  that  the 
central  limit  theorem  is  valid  for  both  scattering  mechanisms,  surface  rough- 
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ness  and  the  bubble  layer,  for  the  geometries  and  all  but  the  highest  wind 
speed  studied  here. 

•  Mean  bistatic  surface  scattering  strength  wcis  calculated  using  the  sonar  equa¬ 
tion  and  compared  to  a  prediction  by  S.T.  McDaniel’s  theory  for  bistatic  sur¬ 
face  scattering  strengths  (15).  The  portion  of  the  theory  for  predicting  scatter¬ 
ing  strengths  dominated  by  bubble  scattering  was  validated  by  the  data  for  a 
number  of  geometries.  This  is  the  first  validation  of  that  portion  of  the  theory. 

•  A  technique  was  developed  for  combining  ocean  surface  scatter  data  for  differ¬ 
ent  grazing  angles  which  increased  the  number  of  data  points  in  the  ensemble. 
With  this  technique,  data  which  were  the  result  of  the  same  dominant  scatter¬ 
ing  mechanism,  the  surface  itself  or  the  sub-surface  bubble  layer,  were  com¬ 
bined  to  estimate  the  statistical  characteristics  of  the  scattering  mechanism. 

1.3  Thesis  Outline 

This  thesis  is  organized  as  follows.  Chapter  2  presents  a  literature  review  of 
the  experimental  work  with  regard  to  both  monostatic  and  bistatic  surface  rever¬ 
beration.  Chapter  3  provides  an  overview  of  the  FLIP  experiment  and  describes  the 
process  by  which  the  data  were  validated  as  homogeneous.  Chapter  4  describes  the 
process  by  which  the  data  were  tested  for  Rayleigh  behavior.  The  results  showed 
the  bistatic  surface  scatter  to  be  Rayleigh  except  at  the  highest  wind  speed  of  7.2 
m/s.  Chapter  5  details  the  calculation  of  the  bistatic  surface  scattering  strength, 
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S S {limine, 'P scat,  4>),  from  the  data  using  the  sonar  equation  for  bistatic  surface  re¬ 
verberation.  Chapter  6  presents  an  overview  of  S.T.  McDamiel's  bistatic  surface 
scattering  strength  theory  and  shows  a  comparison  between  the  theory  and  the 
data.  The  comparison  showed  good  agreement  with  an  rms  deviation  of  about 
3  dB.  Chapter  7  presents  recommendations  for  future  research  in  this  area.  The 
original  contributions  of  this  thesis  are  presented  in  Chapters  4,  5  and  6. 


Chapter  2 


Historical  Background  of  Boundary  Reverberation 

2.1  Experimental  Work  in  Monostatic  Reverberation 

A  number  of  papers  from  the  middle  1950’s  to  the  early  1970’s  provided 
a  great  deal  of  the  knowledge  about  surface  reverberation.  In  general,  they  char- 
acterized  surface  scattering  strength  as  a  function  of  one  or  more  measured  envi* 
ronmental  variables,  such  as  wind  speed  or  wave  height.  Results  were  presented  as 
either  a  set  of  curves  or  an  equation.  However,  variations  of  5  to  10  dB  between 
different  researcher’s  results  were  common.  Nevertheless,  the  data  collected  by 
these  researchers  is  still  valuable  in  identifying  trends  and  general  characteristics 
in  monostatic  surface  reverberation. 

2.1.1  Urick  ie.  Hoover:  1956 

An  experiment  carried  out  by  Urick  and  Hoover  in  1955  [4]  measured  mono¬ 
static  surface  reverberation  for  grazing  angles  between  5®  and  90®,  wind  speeds 
between  3  and  18  knots  and  a  transmit  frequency  of  60  KHz.  The  procedure  for 
extracting  results  from  the  data  was  to  calculate  the  average  scattering  strength 
from  individual  pings  with  the  same  wind  speed  and  grazing  angle.  A  set  of  curves 
were  then  generated  from  the  data,  with  each  curve  corresponding  to  a  particular 


wind  speed,  and  plotted  with  the  grazing  angle  as  the  independent  variable  and 
the  scattering  strength  as  the  dependent  variable.  The  general  trends  shown  in  the 
data  were: 

•  at  low  grazing  angles,  scattering  strength  increases  with  wind  speed, 

•  at  low  grazing  angles  and  higher  wind  speeds,  scattering  strength  becomes  less 
dependent  upon  grazing  angle, 

•  at  grazing  angles  near  normal  incidence,  scattering  strength  decreases  with 
increasing  wind  speed. 

These  trends  are  seen  in  Figure  2.1,  which  shows  monoatatic  surface  scattering 
strength  as  a  function  of  grazing  angle  and  wind  speed.  Figure  2.1  was  generated 
using  the  model  developed  by  S.T.  McDaniel  [15].  More  information  on  this  model 
is  included  in  Chapter  6. 

In  addition,  Urick  and  Hoover  put  forward  several  hypotheses  about  the 
mechamisms  involved  in  sea  surface  scattering.  At  low  grazing  angles,  the  authors 
suggested  that  a  sub-surface  bubble  layer  rather  than  sea  surface  roughness  is 
the  dominant  source  of  scattering.  That  the  scattering  strengths  were  roughly 
independent  of  grazing  angle  is  consistent  with  isotropic  scatterer  by  a  bubble 
layer. 

At  grazing  angles  near  normal  incidence,  the  dominant  echanism  was 
specular  reflection  from  the  sea  surface.  They  theorized  that  the  sea  surface  is 
made  up  of  numerous  facets  whose  orientation  and  distribution  are  dependent  upon 
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wind  speed.  At  very  low  wind  speeds,  the  sea  surface  is  mirror-like,  causing  high 
scattering  strengths  near  normal  incidence  and  low  scattering  strengths  away  from 
normal  incidence.  At  higher  wind  speeds,  the  sea  surface  is  disturbed  and  has  many 
facets  which  scatter  energy  over  a  larger  solid  angle.  The  result  is  lower  scattering 
strengths  near  normal  incidence  and  higher  scattering  strengths  away  from  normal 
incidence. 

The  major  contribution  of  the  work  by  Urick  and  Hoover  was  the  suggestion 
that  the  dominant  mechanism  in  surfau:e  backscatter  depends  upon  the  graizing 
angle  in  the  following  way: 

•  at  low  grazing  angles,  scattering  is  dominated  by  a  sub-surface  scattering  layer, 

•  at  intermediate  grazing  angles,  scattering  is  dominated  by  rough-surface  scat¬ 
tering, 

•  at  high  grazing  angles,  scattering  is  dominated  by  specular  reflection. 

With  the  postulation  that  these  three  mechanisms  were  dominant  at  different  graz¬ 
ing  angles,  other  experiments  sought  to  confirm  these  results. 

2.1.2  Garrison,  Murphy  &c  Potter:  1960 

This  paper  utilized  measurements  taken  by  APL-UW  at  a  transmit  fre¬ 
quency  of  60  KHz  in  Dabob  Bay  during  1956  and  1957  [5]  to  study  the  relationship 
between  surface  reverberation  and  measured  environmental  parameters.  During 
the  experiment,  both  wind  speed  and  wave  height  were  recorded  in  an  effort  to 
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discover  which  of  the  two  was  better  correlated  with  reverberation  levels.  An  in¬ 
teresting  note  about  Dabob  Bay  is  that  the  fetch  from  the  southwest  is  roughly  10 
times  the  length  of  the  fetch  from  the  north.  This  results  in  dramatically  different 
waveheights  for  the  same  wind  speed,  depending  upon  the  wind  direction,  and  was 
particularly  useful  in  isolating  the  effects  of  wind  speed  and  wave  height. 

Transducers  with  circularly  symmetric  beam  patterns  were  deployed  with 
the  axis  of  symmetry  pointed  perpendicular  to  the  surface.  With  this  arrangement, 
the  transmitted  pulse  would  reach  the  surface  and  spread  out  in  an  annular  ring. 
Therefore,  the  reverberation  return  at  any  time  haa  a  corresponding  graizing  angle 
that  ran  be  calculated  from  the  depth  of  the  transducer  and  the  time  after  trans¬ 
mission  of  the  pulse.  The  procedure  in  the  experiment  was  to  transmit  pulses  at  an 
interval  long  enough  to  allow  the  reverberation  from  the  previous  pulse  to  decrease 
to  at  or  below  the  background  noise  level.  The  reverberation  return  was  recorded 
at  a  fixed  time  after  the  pulse  was  transmitted,  corresponding  to  the  grazing  angle 
of  interest.  At  the  same  time,  information  on  the  surface  conditions  were  recorded. 
The  results  of  the  experiment  showed  that  surface  scattering  strength  had  greater 
correlation  with  wind  speed  than  wave  height  at  a  grazing  angle  of  7®.  The  surface 
scattering  strengths  were  more  tightly  grouped  together  at  a  particular  wind  speed 
than  at  a  particular  wave  height. 

Paralleling  the  work  done  by  Urick  and  Hoover  [4],  a  set  of  curves  were 
generated,  each  corresponding  to  a  particular  wind  speed,  with  the  grazing  angle 
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as  the  independent  variable  and  the  surface  scattering  strength  as  the  dependent 
variable.  The  set  of  curves  from  the  APL-UW  data  agree  with  the  curves  from 
Urick  and  Hoover  except  at  low  grazing  angles.  The  APL-UW  curves  show  the 
surface  scattering  strength  falling  off  at  grazing  angles  below  20°  while  the  Urick 
and  Hoover’s  curves  show  the  surface  scattering  strength  remaining  constant  below 
roughly  35°.  This  difference  in  results  may  be  a  consequence  of  the  APL-UW 
data  being  screened  to  omit  points  which  appeared  to  be  contaminated  by  volume 
reverberation,  most  likely  sub-surface  bubble  scatter,  which  were  included  in  the 
Urick  and  Hoover  results. 

2.1.3  Plemmons,  Shooter  &c  Middleton:  1972 

The  work  done  by  Plemmons,  Shooter,  and  Middleton  [7]  at  the  Applied 
Research  Laboratory,  University  of  Texas,  (ARL-UT)  presents  a  statistical  analysis 
of  fresh  water  lake  surface  reverberation.  This  summary  of  their  work  focuses  on 
the  method  by  which  they  tested  the  surface  reverberation  data  for  homogeneity. 

In  all  the  data  sets  collected,  150  pulses  were  transmitted  at  the  operating 
frequency  of  110  KHz  with  an  interval  between  the  pulses  long  enough  for  the 
reverberation  level  to  decay  to  at  or  below  the  background  level.  All  of  the  data 
for  a  particular  pulse  length  and  type  were  collected  over  a  short  period  of  time. 

Quadrature  components  of  the  reverberation  signal  were  grouped  together 
in  ensembles,  each  corresponding  to  the  same  time  after  the  transmission  of  the 
signal.  Each  ensemble  contained  150  values  from  a  distinct  pulse  type  and  delay 


after  transmit.  With  these  ensembles,  estimates  of  certain  characteristics  of  the  re¬ 
verberation  were  calculated,  such  as  the  mean  value  of  the  amplitude  at  a  particular 
time  ti 

1  ^ 

(^(^1  ))enjem6ie  average  ~  "7^  ^i(^l )•  (2.1) 

Before  such  estimates  of  statistical  quantities  were  made,  the  ensembles  were  val¬ 
idated.  Validating  an  ensemble  entails  showing  that  the  ensemble  comes  from  a 
statistically  proper  parent  distribution.  This  involves  testing  the  ensemble  for  ho¬ 
mogeneity. 

The  ensembles  were  tested  for  homogeneity  using  the  Kolmogorov-Smirnov 
(K-S)  two-sample  test,  described  in  detail  by  Siegel  [9),  which  determines  whether 
two  independent  ensembles  have  the  same  parent  distribution.  The  rationale  for 
this  test  is  that  if  two  ensembles,  usually  the  first  and  second  half  of  a  large  en¬ 
semble,  come  from  the  same  parent  distribution,  then  both  will  have  the  same 
cumulative  distribution  function.  The  procedure  is  to  calculate  the  cumulative 
distribution  functions,  cd/’s,  for  each  ensemble  and  find  the  greatest  deviation  be¬ 
tween  the  them.  This  maximum  deviation  is  called  Z^mpie  And  is  compared  with 
Zo,  which  is  determined  from  the  confidence  level,  o,  and  from  the  number  of  sam¬ 
ples  in  the  ensembles.  If  Z,ampie  <  Zo,  the  two  ensembles  are  taken  to  be  from 
the  same  parent  distribution.  If  Z^mpie  >  Zo,  the  ensembles  are  from  different 
parent  distributions.  However,  in  using  this  test,  the  results  of  the  comparison  of 
Zeampie  with  Zo  Are  incorrect  a  x  100%  of  the  time.  This  becomes  important  when 
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large  numbers  of  ensembles  are  tested,  because  some  ensembles  will  fail  the  test  for 
homogeneity  even  when  the  ensembles  are  homogeneous. 

This  paper  presented  a  systematic  method  for  checking  the  validity  of  data 
so  that  results  derived  from  the  data  could  also  be  considered  valid.  This  is  a  key 
step  forward  from  e2U'lier  experiments  which  used  time  averages  over  part  of  a  single 
return  or  averages  over  multiple  returns  without  testing  the  validity  of  the  data. 
This  type  of  analysis  is  used  on  the  FLIP  data  in  Section  3.6. 

2.1.4  FVazer:  1978 

The  paper  by  Marshall  Frazer  [8],  also  from  ARL-UT,  investigated  the  sta¬ 
tistical  properties  of  surface  reverberation  from  a  lake  surface.  The  surface  rever¬ 
beration  data  were  first  tested  for  statistical  validity  and  then  tested  for  normality. 

The  experiment  utilized  a  transducer  located  at  a  depth  of  3.7  meters  with 
the  axis  of  the  main  lobe  intersecting  the  surface  at  a  grazing  angle  of  9".  The 
wind  speed  was  approximately  6.0  m/s  with  a  fetch  of  2.4  km.  A  100-|iscc,  80 
KHz  CW  pulse  was  used  resulting  in  an  scattering  patch  area  of  less  than  1  square 
meter.  This  value  is  important  because  the  assumption  of  Gaussian  characteristics 
in  surface  reverberation  are  based  upon  a  large  number  of  independently  placed 
scatterers.  If  the  reverberation  is  still  Gaussian  with  such  a  small  surface  area 
being  ensonified,  applications  in  which  a  much  larger  surface  area  is  ensonified 
should  also  be  Gaussian. 

The  statistical  validity  of  the  data  was  checked  using  the  one-sample  runs 
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test  and  the  K-S  two-sample  test  as  done  by  Plemmons,  et  al.[7].  With  the  statisti¬ 
cal  validity  of  the  sample  ensembles  established,  the  data  were  tested  for  normality 
using  the  K-S  one-sample  test.  The  mean  of  the  ensemble,  was  calculated 
using  am  ensemble  average,  which  is  an  average  over  all  the  measurements  at  a 
particular  instant  in  time,  t,  after  the  pulse  transmission.  The  variance,  (7^(t),  was 
calculated  in  a  similsu:  manner.  The  mathematical  expressions  for  calculating  the 
ensemble  mean  and  variance  are 

(2.2) 

«=i 

and 

i\t)  =  (2-3) 

With  values  for  fi{t)  and  a  theoretical  Gaussian  (normal)  cumulative  distri¬ 

bution  can  be  calculated  using  the  usual  expression  for  Gaussian  cdf's  with  the 
ensemble  mean  and  variance  introduced 

hhei>rtticai(x{t),V)  =  .  ^  /  exp[-(x  - /i(0)V2^*(0]dx.  (2.4) 

y2jra*(t) 

This  theoretical  cdf  is  then  compared  with  the  cdf  calculated  from  the  ensemble, 
and  the  greatest  deviation  between  the  two  is  the  sample  test  statistic,  Z«amp<e(0* 
The  threshold  value,  Z*,  is  determined  in  the  same  way  as  for  the  K-S  two-sample 
test.  If  Z,amp/e(f)  >  the  ensemble  fails  the  test  and  is  determined  to  be  non- 
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Gaussian  at  the  significance  level  a.  Also,  a  x  100%  of  the  tests  should  fail  even 
if  all  of  the  ensembles  are  Gaussian,  so  judgement  is  called  for  in  analyzing  the 
results  of  the  test. 

The  results  of  the  experiment  indicated  that  the  surface  reverberation  was 
Gaussian  for  the  a  =  0.05  level  of  significance.  However,  it  should  be  noted  that 
some  data  that  contained  both  surface  and  bottom  returns  were  not  Gaussian. 
This  means  that  additional  work  is  needed  to  determine  the  environment  in  which 
reverberation  is  Gaussian  and  more  importantly,  the  environment  in  which  the 
reverberation  is  not  Gaussian.  This  method  for  testing  data  for  normality  was 
used  on  the  FLIP  data  in  Section  4.2. 

2.2  Experimental  Work  in  Bistatic  Reverberation 

Little  work  has  been  done  in  bistatic  scattering  from  the  sea  surface,  and 
the  work  that  has  been  done  covers  a  wide  range  of  frequencies.  Experimental  work 
in  bistatic  surface  scattering  began  in  the  early  1970’s  and  has  continued  to  date, 
but  unlike  the  work  in  monostatic  surface,  scattering  the  mechanisms  responsible 
in  bistatic  surface  scattering  are  less  understood  and  the  data  are  sparse. 

2.2.1  Zornig:  1978 

The  experimental  work  done  by  Zornig  [1]  represents  a  systematic  stud;'  of 
bistatic  surface  scatter.  Working  in  a  tank,  the  wind  conditions  over  the  surface 
were  accurately  controlled  and  repeatable  so  stationary  surface  conditions  could  be 
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maintained  during  measurements.  Operating  at  very  short  wavelengths,  about  1 
mm,  the  Kirchhoff  approximation,  which  is  the  basis  of  much  analytical  work,  could 
be  used.  With  such  a  controlled  environment,  the  effects  of  a  single  parameter  such 
as  wind  direction  could  be  determined.  However,  the  results  of  this  experiment  may 
be  unique  to  the  very  high  frequency,  1.3  MHz,  and  controlled  surface  conditions. 

The  measurements  showed  scattering  strength  to  be  strongly  dependent  on 
wind  direction.  Scattering  strength  was  found  to  be  much  higher  in  the  upwind 
direction  than  in  either  the  downwind  or  crosswind  case.  Other  conclusions  from 
the  results  were  that  scattering  strength  was  not  strongly  dependent  on  either  0,^ 
or  between  12*  and  17®.  This  agrees  with  past  experiments  which  showed 
scattering  strength  to  be  independent  of  grazing  angle  in  the  same  range  of  values. 
Varying  the  acoustic  wavelength  by  changing  the  operating  frequency  from  1.3 
MHz  to  1.1  MHz  produced  little  change  in  the  scattering  strength,  but  variations 
in  surfcice  roughness  produced  large  changes  in  the  scattering  strength. 

Zornig  hypothesized  that  bistatic  surface  scattering  strength  cannot  be 
characterized  by  a  simple  parameter  such  as  the  Rayleigh  roughness  parameter, 
but  may  be  a  function  of  such  descriptors  as  rms  slope  of  the  surface,  probability 
density  of  slopes,  or  two-dimensional  correlation  of  the  surface.  The  work  done  by 
2^rnig  was  a  vaiuaule  step  in  characterizing  the  parameters  that  describe  bistatic 


surface  reverberation. 


2.2.2  Culver  &i  McDaniel:  1991 


The  work  done  by  Culver  and  McDaniel  [3]  incorporated  the  capability  to 
predict  time-varying  reverberation  spectra  into  the  Generic  Sonar  Model  (GSM) 
and  compared  the  predictions  with  bistatic  measurements.  The  new  routine  utilizes 
the  McDaniel  theory  [15]  to  compute  surface  scattering  strength  and  keeps  track  of 
Doppler  shifts  due  to  motion  of  the  source,  receiver,  and  scatterers.  The  McDaniel 
theory  is  described  in  Chapter  6.  Briefly,  it  uses  the  Kirchhoff  approximation  near 
specular  and  the  Rayleigh-Rice  approximation  away  from  specular.  Also,  a  sub¬ 
surface  bubble  layer  incorporated  into  the  model  provides  a  noticeable  effect  at  low 
grazing  angles  and  when  the  number  of  bubbles  per  unit  volume  is  large.  The  effect 
of  the  bubble  layer  is  to  put  a  minimum  value  on  the  surface  scattering  strength 
away  from  the  specular  direction  since  scattering  from  a  bubble  layer  is  isotropic. 

The  bistatic  data  used  for  comparison  were  taken  in  Puget  Sound  with  a 
transmit  frequency  range  of  10  to  30  KHz  and  a  pulse  length  of  0.3  msec.  The  results 
showed  multi-ping  averages  in  a  2.8  KHz  band  around  the  transmit  frequency. 
Overlaid  on  the  experimental  results  were  the  predictions  from  GSM.  The  results 
showed  good  agreement  in  the  portion  of  the  return  dominated  by  the  main  lobe, 
but  are  up  to  30  dB  below  the  measured  values  in  the  portions  of  the  return 
dominated  by  sidelobes.  A  possible  cause  given  for  this  disagreement  is  that  the 
transmit  center  frequency  beam  pattern  was  used  for  the  simulation,  yet  over  the 
bandwidth  of  the  pulse,  the  beam  pattern  may  vary  considerably.  Aside  from  this 
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disagreement,  the  results  show  that  the  surface  return  is  composed  of  two  distinct 
regions.  The  first  region  is  dominated  by  coherent  reflection,  and  is  well  modeled 
by  the  KirchhofF  approximation.  The  second  region,  commencing  roughly  8  to  10 
msec  after  the  pulse  has  arrived,  is  dominated  by  scattering  from  the  rough  surface, 
and  can  be  modeled  by  the  Rayleigh-Rice  approximation.  There  were  no  visible 
effects  of  a  sub-surface  bubble  layer  because  the  data  were  taken  at  large  grazing 
angles  and  low  wind  speeds. 


Chapter  3 


FLIP  Experiment 


3.1  Introduction 

This  chapter  describes  bistatic  surface  scattering  measurements  made  dur¬ 
ing  January  of  1992  from  the  FLIP,  FLoating  Instrument  Platform,  freely  floating 
about  400  nmi  off  the  coast  of  California.  The  experiment  was  planned  and  con¬ 
ducted  by  Peter  Dahl  and  Andrew  Jessup  of  APL-UW,  and  was  primarily  directed 
toward  investigating  the  relationship  between  ocean  surface  conditions,  the  near- 
surface  bubble  layer,  and  high  frequency  acoustic  forward  and  backscatter.  In 
this  chapter,  the  means  by  which  data  were  tested  and  validated  as  homogeneous 
is  described.  The  data's  statistical  characteristics  and  bistatic  surface  scattering 
strengths  will  be  investigated  in  Chapters  4  and  5  respectively. 

3.2  Geometry 

The  geometry  for  the  FLIP  experiment  is  shown  in  Figure  3.1,  which  is 
taken  from  [18].  A  detailed  description  of  FLIP  can  be  found  in  (19). 
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Figure  3.1:  Experimental  Set-up  for  FLIP  Data. 


3.3  Transmitter  and  Receiving  Array 

Three  spherical  sources  were  suspended  below  a  spar  buoy  allowed  to  drift 
approximately  1000  meters  from  FLIP.  The  sources  were  at  depths  of  27,  57  and 
147  m.  Any  one  of  the  three  sources  could  be  selected  to  transmit.  The  transmitted 
pulse  was  received  by  horizontal  and  vertical  line  arrays  mounted  on  FLIP’S  hull  at 
a  depth  of  66.3  m.  Each  line  array  was  composed  of  84  omnidirectional  elements 
with  0.0153  m  spacing.  The  84  omnidirectional  elements  were  grouped  into  4  staves, 
each  with  21  elements.  The  signals  from  each  stave,  eight  in  all,  were  shifted  down 
to  a  center  frequency  of  5  KHz  and  sampled  at  20  KHz. 

The  water  depth  at  the  experiment  site  was  approximately  4000  meters. 


The  depth  of  the  water  kept  the  returns  at  the  line  arrays  limited  to  just  the  direct 
path,  the  specular  reflection  from  the  surface  and  the  surface  scatter. 

3.3.1  Processing  the  Line  Array  Data 

The  objective  in  processing  or  combining  line  array  stave  data  was  to  re¬ 
solve  a  single,  small  scattering  patch  on  the  surface  so  that  the  bistatic  surface 
scattering  strengths  derived  from  the  data  could  be  related  to  particular  geometri¬ 
cal  parameters.  This  section  discusses  four  possible  methods  of  processing  the  line 
array  data  and  shows  that  using  the  horizontal  array  alone  proved  to  be  best  choice 
for  resolving  the  scattering  patch. 

The  four  possible  methods  of  processing  the  line  array  data  were  to  form 
a  Mills  Cross  array  of  the  horizontal  and  vertical  arrays,  form  a  sum  array  of 
the  horizontal  and  vertical  arrays,  use  the  vertical  array  alone  or  use  the  horizontal 
array  alone.  The  four  resulting  beam  pattern  plots  are  shown  in  Figures  3.2  through 
3.5. 

Summing  the  returns  from  the  four  vertical  staves  and  multiplying  by  the 
sum  of  the  four  horizontal  staves  outputs  produces  a  Mills  Cross  array  [17].  Pro¬ 
cessing  the  output  of  the  line  arrays  in  this  manner  produces  a  beam  pattern  that 
is  the  product  of  the  beam  patterns  of  the  horizontal  and  vertical  line  arrays.  A 
plot  of  the  Mills  Cross  array  beam  pattern  is  included  in  Figure  3.2.  There  is  a  0 
dB  peak  at  0"  bearing  and  0"  elevation  and  high  sidelobes  in  the  horizontal  and 
vertical  planes. 
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Figure  3.2:  Mills  Cross  Array  Beam  Pattern. 

Adding  all  eight  staves  produces  a  sum  array  with  a  beam  pattern  con¬ 
taining  very  high  sidelobes  in  the  horizontal  and  vertical  planes.  The  sum  array 
improves  the  signal-to-noise  ratio  when  compared  to  either  of  the  line  arrays  treated 
separately,  in  contrast  to  the  Mills  Cross  array,  which  was  found  to  suffer  from  a 
lower  signal-to-noise  ratio  when  compared  to  either  of  the  line  arrays  treated  sep¬ 
arately.  A  plot  of  the  sum  array  beam  pattern  is  included  in  Figure  3.3  and  shows 
the  same  general  shape  as  the  Mills  Cross  array  but  with  higher  sidelobes. 

Summing  the  staves  separately  for  each  line  array  produces  the  expected 
beam  pattern.  In  Figure  3.4,  the  vertical  line  array  beam  pattern  is  shown.  Note 
that  there  is  a  plane  of  ambiguity  at  0”  elevation.  In  Figure  3.5  the  horizontal  line 
array  beam  pattern  is  shown  with  the  ambiguity  plane  at  0*  bearing.  The  ridges  at 
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Figure  3.3:  Sum  Array  Beam  Pattern. 


±90*  elevation  are  actually  points,  in  that  bearing  angle  has  no  meaning  at  these 
elevation  angles. 

The  detrimental  effect  of  the  ambiguity  plane  at  0*  elevation  in  ail  the  beam 
patterns  except  that  of  the  horizontal  array  alone  proved  to  be  the  reason  for  using 
the  horizontal  array  alone. 

A  graphical  representation  of  this  problem  is  given  next.  In  Figure  3.6, 
the  scattering  patches  are  shown  at  675,  700,  725,  750,  775  and  800  msec  after 
transmit  for  the  horizontal  array  alone.  Here  the  scattering  patch  is  defined  as 
all  those  points  within  3  dB  of  the  maximum  of  the  beam  pattern  gains  minus 
the  transmission  and  attenuation  losses.  In  Figure  3.6,  the  source  is  located  at 
(•500.0,  0.0)  and  the  receiver  at  (500.0,  0.0).  The  scattering  patch  for  time  675 


Figure  3.4:  Vertical  Line  Array  Beam  Pattern. 
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Figure  3.6:  Scattering  Patches  (within  3  dB)  for  Horizontal  Array  at  Times  675  to 
800  msec  in  25  msec  intervals. 

msec  is  located  at  approximately  (350.0,  0.0).  The  scattering  patches  for  times  700 
to  800  msec  start  close  to  the  receiver  and  move  away  from  the  receiver  for  longer 
times.  Figure  3.7  is  analogous  to  Figure  3.6  except  that  the  scattering  patches  are 
defined  as  all  those  points  within  15  dB  of  the  maximum  value.  The  increase  in 
the  scattering  patch  area  between  Figure  3.7  and  Figure  3.6  is  confined  to  the  area 
around  the  same  physical  location  on  the  surface  at  each  time  and  is  up  to  a  factor 
of  three  times  larger  except  at  time  675  msec. 

In  contrast,  Figure  3.8  shows  the  scattering  patch  defined  as  those  points 
within  3  dB  of  the  maximum  for  the  same  times  as  before  but  for  the  sum  beam 
pattern.  The  scattering  patch  at  three  of  the  six  times  is  not  contiguous.  In 
Figure  3.9,  the  sum  beam  15  dB  scattering  patches  are  shown.  For  all  times,  the 
scattering  patches  are  not  contiguous  and  it  is  impossible  to  relate  a  single  geometry 
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Figure  3.9:  Scattering  Patches  (within  15  dB)  for  Sum  Array  at  Times  675  to  800 
msec  in  25  msec  intervals. 

to  a  bistatic  surface  scattering  strength  measurement.  This  is  also  true  with  all 
array  configurations  which  include  the  vertical  array  because  of  its  ambiguity  plane 
at  0"  elevation. 

3.4  Environmental  Parameters  Measured 

Environmental  parameters  measured  during  the  FLIP  experiment  include 
wind  speed,  wind  direction,  and  air  temperature.  The  sound  speed  profile  was 
measured  twice  daily  during  the  entire  experiment.  The  environmental  data  is 
used  to  investigate  the  effects  of  the  environment  on  bistatic  scattering  strengths. 
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Table  3.1:  Transmitter  Information  for  Bistatic  Surface  Scattering  Portion  of  the 
FLIP  Experiment. 


Run 

Freq  (KHz) 

Pulse  (ms) 

SL  (dB) 

#  of  Pings 

fwd-las.  102 

30 

8 

188.2 

50 

fwd-las.  103 

30 

8 

188.2 

fwd-las.  104 

30 

4 

189 

25 

fwd-las.  105 

30 

4 

189 

25 

fwd-las.  106 

30 

4 

189 

25 

fwd-las.  107 

30 

4 

25 

fwd-las.  108 

40 

10 

193 

25 

fwd-las.  112 

30 

1 

198 

fwd-las.  112 

40 

1 

197 

20 

fwd-las. 116 

30 

1 

195 

25 

fwd-las.  117 

30 

1 

195 

25 

3.5  Acoustic  Data 

The  bistatic  surface  scattering  measurements  made  during  the  FLIP  exper¬ 
iment  are  described  in  the  Tables  3.1,  3.2  and  3.3.  Run  refers  to  the  type  of  data, 
where  fwd-las  means  forward  line  array  sonar,  and  the  number  indicates  a  particu¬ 
lar  data  set.  The  elevation  and  bearing  of  the  line  arrays  is  given  in  Table  3.3  and 
is  expressed  in  degrees.  A  bearing  angle  of  0"  corresponds  to  the  direction  of  the 
source  and  a  elevation  angle  of  0°  corresponds  to  horizontal.  The  time  delay  Td, 
is  the  record  start  time  referenced  to  the  transmission  of  the  pulse,  and  the  time 
sample  T,  is  the  length  of  the  record  time.  Both  are  expressed  in  msec.  The  wind 
speed  during  each  run  is  given  in  Table  3.4.  At  the  experiment  site,  the  received 
signals  were  shifted  down  to  a  center  frequency  of  5  KHz  with  a  bandwidth  of  10 


Table  3.2:  Transmitter  Location  during  Bistatic  Surface  Scattering  Runs. 


Run 

Source  Depth  (m) 

Source  Range  (m) 

fwd-las.l02 

147 

977 

fwd-las.103 

147 

977 

fwd-las.104 

57 

999 

fwd-las.105 

57 

1002 

fwd-las.108 

147 

714 

fwd-las.112 

147 

980 

fwd*Ias.ll6 

147 

715 

fwd-las.117 

147 

712 

Table  3.3:  Receiver  Information  for  Bistatic  Surface  Scattering  Portion  of  FLIP 
Experiment. 


Run 

Rec.  Elev. 

Rec.  Brng. 

mm 

mm 

Gain  (dB)  (| 

fwd*las.l02 

10 

55 

600 

200 

fwd-las.103 

20 

55 

600 

200 

50 

fwd-las.104 

15 

83 

600 

200 

50 

fwd-las.105 

5 

83 

600 

200 

50 

fwd-las.108 

15 

8 

450 

100 

38 

fwd-las.112  (30KHz) 

2.5 

55 

650 

100 

50 

fwd-las.112  (40KHz) 

2.5 

55 

650 

100 

50 

fwd-b’^.llS 

2.5 

30 

475 

100 

56 

fwd-las.117 

2.5 

60 

475 

100 

56 

30 


Table  3.4:  Measured  Wind  Speeds  for  Bistatic  Surface  Scattering  Measurements. 


Runs 

Date 

Time 

Wind  Speed  (m/s) 

fwd-las.l02 

1/17/92 

1535 

7.2 

fwd-las.103 

1/17/92 

1536 

7.2 

fwd-las.104 

1/18/92 

1930 

5.5 

fwd-l2is.l05 

1/18/92 

2036 

5.8 

fwd-las.108 

1/21/92 

1208 

1.4 

fwd-las.ir2 

1/22/92 

2005 

5.3 

fwd-las.ll6 

1/24/92 

1140 

5.3 

fwd-las.117 

1/24/92 

1154 

4.9 

KHz  and  then  sampled  at  20  KHz.  The  data  were  stored  in  two  byte  binary  format, 
time  multiplexed  for  the  8  staves. 

3.6  Testing  for  Homogeneity 

Before  investigating  the  statistics  of  the  time  series,  the  stationarity  of  the 
data  was  verified  using  the  two-sample  Kolmogorov-Smirnov  test.  The  procedure 
used  was  the  same  as  detailed  in  Section  2.1.3.  Ageun,  the  purpose  of  this  test  is 
to  determine  whether  the  statistical  characteristics  of  the  reverberation  changed 
between  the  first  half  and  second  half  of  the  transmitted  pulses.  This  is  an  impor¬ 
tant  step  because  the  test  for  normality  is  only  valid  if  the  statistics  of  the  entire 
ensemble  are  homogeneous. 

The  procedure  for  checking  the  data  for  homogeneity  was  to  divide  the 
data  into  a  first  half  and  a  second  half  and  then  calculate  a  cumulative  distribution 


Table  3.5:  Test  Thresholds  for  K-S  Test  for  Homogeneity. 


Run 

mm 

ni 

nj 

fwd-las.102 

.345 

25 

25 

fwd-las.l03 

.345 

25 

25 

fwd-las.104 

.488 

12 

13 

fwd-las.105 

.488 

12 

13 

fwd-las.108 

.488 

12 

13 

fwd-las.112 

.545 

10 

10 

fwd-las.116 

.488 

12 

13 

fwd-las.117 

.488 

12 

13 

function,  cd/,  at  each  time  for  each  half.  Then,  one  finds  the  maximum  deviation 
between  the  two  cdf’s  at  each  time  and  uses  this  value,  Ztampte^  as  the  test  sample 
statistic.  This  value,  Zaampie,  is  compared  with  the  test  threshold,  Zo,.  The  test 
threshold  was  calculated  using  the  formula 


Za  =  1.22, 


fni  +  n2 

ninj 


where  ni  and  are  number  of  data  points  in  the  two  ensembles  being  compared 
[9].  In  this  application  of  the  K-S  test,  ni  is  equal  to  the  number  of  pulses  in  the 
first  half  of  the  ensemble,  is  equal  to  the  number  of  pulses  in  the  second  half  of 
the  ensemble  and  the  confidence  level,  q,  is  equal  to  0.1.  The  values  of  Za  used  in 
the  test  for  homogeneity  as  well  as  the  values  of  ni  and  n2  are  shown  in  Table  3.5. 
Test  results  are  displayed  in  Figures  3.10  through  3.12.  It  should  be  noted  that 
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when  testing  large  numbers  of  ensembles  it  is  acceptable  for  up  to  (100  x  a)% 
of  the  test  statistics  to  exceed  the  threshold  before  the  data  set  fails  the  test. 
Figures  3.10  through  3.12  show  that  all  of  the  runs  pass  the  test  for  homogeneity 
using  the  Kolmogorov-Smirnov  test  with  a  confidence  level  of  q  =  0.10  in  that  no 
more  than  10%  of  the  test  statistics  exceed  the  threshold  for  any  run. 

3.7  Conclusion 

By  testing  data  for  homogeneity,  one  can  avoid  the  mistake  of  averaging 
over  a  period  of  time  during  which  the  underlying  statistics  change  significantly. 
This  is  an  important  consideration  because  derived  values  such  as  the  bistatic 
surface  scattering  strength  are  inherently  connected  to  statistical  parameters  such 
as  the  variance  of  the  reverberant  time  series.  If  the  v^iance  changes  during  the 
sequence  of  transmitted  pulses,  then  an  estimate  of  the  bistatic  surface  scattering 
strength  averaged  over  the  sequence  will  be  inherently  flawed. 

With  all  the  bistatic  runs  validated  as  homogeneous  at  a  confidence  level  of 
a  =  0.10,  they  could  then  be  tested  for  Rayleigh  behavior  and  the  bistatic  surface 
scattering  strength  calculated  as  discussed  in  Chapters  4  and  5. 
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Figure  3.12:  Testing  for  Homogeneity  •  fwd-las.117.  versus  time  (msec)  is 

plotted.  The  horizontal  line  is  Zq. 


Chapter  4 


FLIP  Data  -  Statistics 


4.1  Introduction 

In  this  chapter,  the  results  of  tests  for  Rayleigh  behavior  using  both  the 
Kolmogorov- Smirnov  (K-S)  two  sample  test  and  a  combinational  technique  de¬ 
signed  to  increase  the  number  of  points  in  the  ensembles  are  described.  The  com¬ 
binational  technique  resolved  detail  in  the  tail  of  the  distribution  and  revealed 
deviations  from  Rayleigh  behavior  at  the  7.2  m/s  wind  speed  which  were  not  seen 
using  the  Kolmogorov-Smirnov  test  for  Gaussian  behavior. 

4.2  K-S  Test  for  Normality 

The  two-sample  K-S  test  can  be  applied  in  a  somewhat  different  manner 
to  that  described  in  Chapter  3  to  determine  whether  the  data  exhibit  Gaussian 
statistical  characteristics.  For  this  application,  the  cdf  of  all  of  the  pulses  at  a 
particular  time  is  compared  with  the  cdf  of  a  Gaussian  random  sequence  of  1000 
points  having  the  same  mean  and  variance  as  that  of  the  pulses.  The  test  statistic 
is  taken  to  be  the  maximum  deviation  between  the  two  cd/’s,  and  is  compared  to 
a  test  threshold  as  before.  Table  4.1  shows  the  values  of  the  test  threshold  (with 
a  =  0.1)  and  the  values  of  ni  and  nj  used  in  the  test  for  normality. 
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Table  4.1;  Test  Thresholds  for  K-S  Test  for  Norm«ility. 


Run 

mm 

ni 

712 

fwd-las.102 

.176 

50 

1000 

fwd-las.103 

.176 

50 

1000 

fwd-las.104 

.247 

25 

1000 

fwd-las.105 

.247 

25 

1000 

fwd-las.108 

.247 

25 

1000 

fwd-ias.112 

.275 

20 

wm\ 

fwd-las.116 

.247 

25 

DUSII 

fwd-las.117 

.247 

25 

1000 

The  results  of  the  test  for  normality  are  shown  in  Figures  4.1,  4.2  and  4.3. 
All  runs  passed  the  test  in  that  few  points  exceed  the  threshold.  However,  there 
are  two  deficiencies  of  this  result.  First,  each  value  of  Ztampit  is  based  upon  at 
most  50  reverberation  data  samples.  This  is  too  few  samples  to  observe  deviations 
in  the  tails  of  the  distribution,  where  events  occur  only  once  in,  say  10^  samples. 
Second,  the  K-S  test  looks  at  the  entire  distribution,  whereas  we  are  most  interested 
in  the  shape  of  the  distribution  tails.  In  the  next  section,  the  data  statistics  are 
investigated  in  a  manner  which  remedies  both  deficiencies. 


4.3  Testing  the  Envelope  for  Rayleigh  Behavior 

In  order  to  increase  the  number  of  samples  in  the  test  for  normality,  data 
from  different  times  were  normalized  and  combined  into  a  single  ensemble.  Por¬ 
tions  of  the  time  series  in  which  a  single  mechanism  was  dominant,  either  surface 
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Test  far  Normal  itu  -  fud-lcs.117 


Figure  4.3:  Testing  for  Normality  -  fwd>las.ll7.  Z^mpte  versus  time  (msec)  is 
plotted.  The  horizontal  line  is  Z^. 

scattering  or  bubble  scattering,  were  combined.  Based  upon  work  by  Urick  and 
Hoover  [4]  and  Culver  and  McDaniel  [3],  the  initial  reverberation  return  is  due  to 
high  angle  scattering  and  is  governed  by  surface  scattering.  At  later  times,  bubble 
scattering  is  dominant.  The  time  at  which  the  crossover  occurs  was  determined  us¬ 
ing  the  McDaniel  theory  for  predicting  bistatic  surface  scattering  strength,  which 
separately  predicts  the  effects  of  surface  roughness  and  the  bubble  layer.  Table  4.2 
shows  the  time  period  and  available  data  points  during  which  each  mechanism  was 
dominant.  A  portion  of  each  run  is  dominated  by  surface  scattering  and  a  por¬ 
tion  is  dominated  by  bubble  scattering,  except  for  run  fwd-las.108  which  is  entirely 
dominated  by  surface  scattering. 

The  equivalency  of  testing  for  Rayleigh  behavior  in  the  magnitude  of  the 
complex  envelope  and  testing  for  normal  behavior  in  the  time  series  is  discussed  in 
[10].  Basically,  if  the  real  and  imaginary  parts  of  a  complex  time  series  are  normally 
distributed  random  variables,  the  magnitude  of  the  complex  envelope  will  have  a 
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Table  4.2;  Receiver  Information  for  FLIP  Data. 


Run 

Begin  Time 

End  Time 

#  points 

Mechanism 

fwd-las.102 

.675 

.715 

250 

Surface 

fwd*las.l02 

.715 

.800 

550 

Bubbles 

fwd-las.103 

.675 

.724 

350 

Surface 

fwd-las.103 

.725 

.800 

500 

Bubbles 

fwd-las.104 

.675 

.724 

325 

Surface 

fwd-las.104 

.725 

.800 

475 

Bubbles 

fwd-Ias.105 

.680 

.722 

275 

Surface 

fwd-las.105 

.723 

.800 

500 

Bubbles 

fwd-las.108 

.510 

.550 

100 

Surface 

fwd-las.112  (30KHz) 

.675 

.721 

920 

Surface 

fwd-las.112  (30KHz) 

.722 

.750 

560 

Bubbles 

fwd-las.112  (40KHz) 

.675 

.721 

920 

Surface 

fwd-las.112  (40KHz) 

.722 

.750 

560 

Bubbles 

fwd-las.116 

.502 

.531 

725 

Surface 

fwd-las.116 

.532 

.575 

1075 

Bubbles 

fwd-las.117 

.502 

.537 

875 

Surface 

fwd-las.117 

.538 

.575 

925 

Bubbles 
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Rayleigh  distribution. 

Once  a  basebanded,  filtered  complex  time  series  was  calculated,  the  linear 
magnitude  was  taken  and  subsampted  by  the  length  of  the  pulse  to  provide  indepen¬ 
dent  sample  of  the  magnitude  of  the  complex  envelope.  The  data  for  each  time  were 
normalized  to  a  mean  of  10  and  a  variance  of  1.  Setting  the  mean  to  10  ensured  all 
the  samples  were  positive.  The  normalized  data  for  each  time  were  then  combined 
with  data  from  other  times  with  the  same  dominant  scattering  mechanism.  The 
Pfa  was  generated  by  histogramming  the  ensemble  and  taking  one  minus  the  cdf. 

The  results  are  plotted  on  a  log-log  scale  with  normalized  envelope  level  on 
the  horizontal  axis  amd  the  probability  of  false  alarm,  Pfa-,  on  the  vertical  aocis.  The 
Pfa  highlights  deviations  from  Rayleigh  behavior  in  the  tails  of  the  distribution. 
An  example  of  this  type  of  analysis  used  on  bottom  backscatter  is  given  in  an  article 
by  Boehme  and  Chotiros  [20]. 

4.3.1  Surface  Roughness  Scattering 

The  Pfa  results  are  presented  in  order  of  decreasing  wind  speed  since  the 
wind  speed  is  the  driving  force  behind  surface  disturbance.  The  order  of  the  runs  is 
given  in  Table  4.3,  which  shows  the  run,  the  wind  speed  and  the  transmit  frequency. 

Note  that  all  of  the  data  presented  are  at  30KHz,  except  for  fwd-las.112,  which 
is  at  30  and  40  KHz,  and  fwd-las.l08,  which  is  at  40  KHz.  Only  the  surface 
scattering  statistics  for  run  fwd-las.102  showed  meaningful  deviations  from  the 
Rayleigh  distribution  based  on  a  visual  inspection  of  the  plots.  Fwd-las.102  shows 
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PFR  for  fwdlos.  103  (tintes  675  to  72d  msec)  -  Surface  ScoTter:nc 


20  lag  of  Envalopt  LmI  (Narnnilzid  to  vorlcnca  1  with  a  ffl«n  of  IS) 

Figure  4.5:  Surface  Scattering  Statistics  fwd>las.l03. 

a  somewhat  higher  probability  of  false  alarm  than  Rayleigh  in  the  tails  of  the 
distribution.  This  is  not  an  unexpected  result  for  a  relatively  high  ^vind  speed.  This 
run  occurred  at  the  highest  wind  speed  of  all  the  runs,  7.2  m/s.  It  is  interesting 
to  note  that  run  fwd-ias.i03,  taken  only  minutes  later  with  the  same  wind  speed 
and  an  elevation  of  20*  instead  of  10*,  showed  no  significant  deviation  from  the 
Rayleigh  distribution. 

4.3.2  Bubble  Li^er  Scattering 

As  in  Section  4.3.1,  the  Pfa  results  for  scattering  from  the  bubble  layer  are 
presented  in  order  of  decreasing  wind  speed.  Significant  deviation  from  the  Rayleigh 
distribution  occur  only  for  runs  fwd'las.102  and  103,  the  runs  that  occurred  at  the 
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PFfl  for  fwdlas.105  (‘times  660  to  722  msec)  -  Surface  Scotte^'irg 


20  log  of  Enwatope  Levtl  (Nomollnd  to  wcrlonco  I  with  a  neon  of  101 

Figure  4.6:  Surface  Scattering  Statistics  fwd*las.l05. 

PFfi  for  fwd\os.l04  (times  675  to  72fl  msec)  -  Surfoce  Scottering 


20  tog  of  Envalop*  Laval  (^knlaliz■d  to  vo-lonca  1  with  a  maon  of  10) 

Figure  4.7:  Surface  Scattering  Statistics  fwd-las.104. 
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PFfi  for  fwdlas.116  (tirnes  502  to  531  msec)  *  Surfoce  Scatterfng 


20  log  of  Envelope  Level  (hkrmllzed  to  variance  1  with  o  moi  of  10) 

Figure  4.8:  Surface  Scattering  Statistics  fwd«Ias.ll6. 

PFfi  for  fwdlas.112  -  30KHz  (times  675  to  721  msec)  •  Sirfoce  Scattering 


Figure  4.9:  Surface  Scattering  Statistics  fwd*las.ll2*30KHz. 
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20  log  of  Enwaiopo  Laval  (Normal  Izad  to  wtrlonca  1  with  o  moon  of  10) 


Figure  4.10:  Surface  Scattering  Statistics  fwd-la8.112-40KH2. 
PFfl  for  fwdlas.117  (limes  502  to  537  msec)  -  Surface  Scottering 


16.0  16.5  17.0  17.5  18.0  18.5  19.0  19.5  20.0  20.5  21.0  21.5  22.0  22.5  23.0  23.5  24.0 
20  log  of  Envalopa  Laval  (Nrmollzod  to  vorlonca  1  with  o  moon  of  10) 


Figure  4.11:  Surface  Scattering  Statistics  fwd-las.117. 


16.0  16.5  17.0  17.5  18.0  18.5  19.0  19.5  20.0  20.5  21.0  21.5  22.0  22.5  23.0  23.5  24.0 
20  log  of  Enwolopo  Uaval  (hkniolizid  to  vorionco  1  with  a  iwcn  of  10) 


Figure  4.12:  Surface  Scattering  Statistics  fwd-las.lOS. 
PFP  for  fwdlas.102  (ilmes  715  to  800  msec)  *  Sufafale  Scattering 


16.0  16.5  17.0  17.5  18.0  18.5  19.0  19.5  20.0  20.5  21.0  21.5  22.0  22.5  23.0  23.5  24.0 
20  log  of  Envelopo  Lawol  (^tanlallz■d  to  vorlonca  1  with  o  naan  of  10) 

Figure  4.13:  Bubble  Scattering  Statistics  {wd-las.l02. 
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PFfl  for  fwdlas.103  (titnes  725  to  800  msec)  -  Bubble  Scattering 


20  log  of  Envelope  Level  (Nomollzed  to  vorlonce  1  with  o  tneon  of  10) 

Figure  4.14:  Bubble  Scattering  Statistics  fwd'las.lOS. 

highest  wind  speed,  7.2  m/s.  The  shape  of  the  probability  of  false  alarm  curves 
for  both  runs  fwd-Us  102  and  fwd*las.l03  show  the  same  trend.  They  are  slightl 
under  and  then  slightly  over  the  Rayleigh  distribution  further  into  the  tail  of  the 
distribution^  The  rest  of  the  plots  for  lower  wind  speeds  show  little  deviation  from 
the  Rayleigh  distribution. 

4.4  Summary  and  Conclusions 

While  all  of  the  data  sets  passed  the  K-S  test  for  normidity,  this  is  not 
considered  conclusive  because  only  25  to  50  points  were  available  for  an  ensemble 
and  because  the  K*S  test  does  not  focus  on  the  tails  of  the  distributions  where 
non-Gaussian  behavior  is  of  most  interest.  In  order  to  increase  ensemble  size,  data 


16.e  16.5  17.0  17.5  18.0  18.5  19.0  19.5  20.0  20.5  21.0  21.5  22.0  22.5  23.0  23.5  24.0 
20  log  of  Envolopo  Lavol  (Nonnltzod  to  vorino  1  with  o  mocr  of  10) 

Figure  4.15:  Bubble  Scattering  Statistics  fwd-las.  105. 


PFfl  for  fvdlos.  104  (times  725  to  800  msec)  -  Bubble  Scotter jng 
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20  log  of  Enwelopt  Levol  Mmollad  to  vorionoo  1  with  a  neon  of  10) 
Figure  4.16:  Bubble  Scattering  Statistics  fwd-las.  104. 
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20  log  of  Enwelope  Level  (Nornatleed  to  variance  1  with  a  neon  of  10) 

Figure  4.17;  Bubble  Scattering  Statistics  fwd'las.llG. 

PFfl  for  fwdlas.112  -  30KHz  (times  722  to  750  msec)  *  Btbble  Scottering 


20  log  of  Envelope  Level  (Narmollzed  to  variance  1  with  a  mean  of  10) 

Figure  4.18:  Bubble  Scattering  Statistics  fwd*la8.112-30KHz. 
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Figure  4.19:  Bubble  Scattering  Statistics  fwd-las.ll2*40KHz. 
PFfi  for  fwdlos.117  (times  538  to  575  msec)  -  Bubble  Scottering 
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20  log  of  Envelopo  Lavol  (Nomallxod  to  worlonca  1  with  a  moon  of  10) 

Figure  4.20:  Bubble  Scattering  Statistics  fwd-las.llT. 


53 


from  different  times  were  normalized  and  combined  for  periods  with  a  common 
dominant  mechanism,  either  rough  surface  or  bubble  layer  scattering,  for  each  data 
run.  A  complex  Gaussian  distributed  process  will  have  a  Rayleigh  envelope,  and 
the  probability  of  false  alarm,  was  used  to  investigate  non-Gaussian  behavior 
in  the  distribution  tails.  It  was  found  that  for  both  mechanisms,  surface  roughness 
and  bubble  scatter,  at  the  highest  wind  speed  of  7.2  m/s,  there  were  meaningful 
deviations  from  Rayleigh  behavior.  The  deviations  are  well  below  the  threshold  of 
the  K-S  test  with  a  confidence  level  of  a  =  0.1,  which  shows  the  limitations  of  the 
K-S  test  in  detecting  deviations  in  the  tails  of  distributions. 


Chapter  5 


FLIP  Data  -  Bistatic  Scattering  Strengths 

5.1  Introduction 

This  chapter  discusses  how  the  FLIP  experiment  data  listed  in  Table  3.1 
were  processed  to  yield  estimates  of  bistatic  surface  scattering  strength  for  various 
geometries.  Only  data  during  the  portion  of  the  time  series  in  which  the  scattering 
patch  area  was  a  single,  relatively  small  area  on  the  surface  could  be  processed.  In 
Chapter  6,  the  scattering  strength  estimates  aure  compared  to  the  McDaniel  model. 

5.2  Finding  the  Scattering  Patch  Area 

A  crucial  first  step  in  calculating  bistatic  surface  scattering  strength  from 
the  FLIP  data  is  to  estimate  the  scattering  patch  size  and  location.  Appendix  A 
provides  the  basics  of  bistatic  scattering  patch  geometry.  A  computationally  inten¬ 
sive  grid  search  method  was  chosen  to  find  the  scattering  patch.  The  alternative 
method  required  solving  a  system  of  non-linear  equations  to  arrive  at  analytic  ex¬ 
pressions  for  the  equi-time  ellipses  on  the  surface,  at  which  point  one  would  still 
have  to  incorporate  the  effects  of  beam  gains,  transmission  losses  and  attenuation 
losses  to  arrive  at  an  estimate  of  the  scattering  patch  area.  The  method  used  in  this 
thesis  required  calculating  the  ensonified  area  at  each  instant  of  time,  then  normal- 
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izing  by  the  transmission  losses,  attenuation  due  to  absorption,  and  beam  patterns. 
The  FORTRAN  code  used  to  estimate  the  bistatic  surface  scattering  strength  is 
included  in  Appendix  B.  The  program  places  the  source  and  receiver  on  the  x  axis 
symmetrically  about  the  origin  and  then  searches  a  grid  of  equally  spaced  points  on 
the  surface  to  find  all  the  points  which  are  within  the  equi-time  ellipses  delimited 
by  the  time  of  interest  and  the  pulse  length.  For  ezwrh  point  within  the  equi-time 
ellipse,  the  transmission  losses,  the  beam  pattern  gains  and  the  attenuation  losses 
are  calculated  and  summed  into  the  variable  test.  The  program  finds  the  point 
which  has  the  highest  value  of  test,  which  corresponds  to  the  point  with  the  lowest 
losses.  All  the  points  whose  value  of  test  are  within  3  dB  of  the  maximum  value 
of  test  are  counted  as  part  of  the  scattering  patch,  and  the  area  of  the  scattering 
patch  is  taken  to  be  the  number  of  points  times  the  grid  spacing  squared.  The  size 
of  the  grid  spacing  is  then  reduced  by  a  factor  of  two  and  the  process  repeated 
to  arrive  at  another  estimate  of  the  scattering  patch  area.  If  the  second  estimate 
satisfies  the  inequality 

1.25  X  Afpacin#/}  ^  Ajpocinj  ^  0*8  X  Afpacing/2  (5.1) 

then  the  area  calculated  using  the  smaller  spacing  is  used.  If  this  inequality  is  not 
satisfied,  the  spacing  is  again  reduced  by  a  factor  of  two  and  the  area  recalculated. 
The  process  stops  when  Equation  5.1  is  satisfied. 

With  this  estimate  of  the  scattering  patch  area.  A,  bistatic  surface  scatter- 


ing  strength  was  calculated  using  the  sonar  equation  solved  for  the  bistatic  surface 
scattering  strength. 

SS  —  RL  +  T L\j^  +  T* —  SL  —  10  log  A  —  BG  +  oiRia  ^  qRza  (^■2) 

The  terminology  used  in  the  sonar  equation  above  is  as  follows. 

•  SS  is  the  bistatic  surface  scattering  strength  (dB). 

•  RL  is  the  received  level  at  the  hydrophone  (dB//^Pa). 

•  TLxa  is  the  transmission  loss  from  the  source  to  the  scattering  patch  (dB). 

•  TL2A  is  the  transmission  loss  from  the  scattering  patch  to  the  receiver  (dB). 

•  SL  is  the  source  level  (dB//^Pa  @  Im). 

•  is  the  ensonihed  area  (m^). 

•  BG  is  the  beam  gain  (dB). 

•  a  is  the  attenuation  coefficient  (dB/m). 

•  Ria  is  the  distance  from  the  source  to  the  scattering  patch  (m). 

•  RiA  is  the  distance  from  the  scattering  patch  to  the  receiver  (m). 

However,  before  this  process  can  be  implemented,  a  value  for  the  received  level, 
RL,  must  be  calculated. 


There  is  one  possible  cause  of  error  that  should  be  mentioned  for  complete¬ 
ness.  The  assumption  was  made  that  an  iso-velocity  sound  profile  was  a  reasonable 
approximation  to  the  medium  even  though  there  were  available  sound  velocity 
profiles  taken  at  the  experiment  site.  To  check  the  effect  of  the  measured  sound 
velocity  profiles  on  propagation  over  the  range  from  the  source  to  the  receiver, 
ray  traces  were  run  on  the  sound  'velocity  profiles  to  look  for  any  large  deviations 
from  a  iso-velocity  profile  approximation.  In  all  of  these  comparisons,  the  amount 
of  refraction  in  the  ray  traces  was  minimal  over  the  range  of  propagation  in  the 
experiment,  less  than  1  km.  This  was  quantified  by  looking  at  the  arrival  angles  of 
the  ray  traces  compared  to  those  using  an  iso-velocity  approximation.  These  arrival 
angles  differed  by  at  most  several  degrees  which  is  comparable  to  the  accuracy  of 
the  measured  elevation  angle  and  bearing  angle  for  the  received  line  arrays. 

5.3  Calculating  the  Received  Level 

With  the  stationarity  of  the  data  verified  in  Section  3.6,  the  data  were 
processed  to  arrive  at  the  received  level.  The  average  spectrum  of  a  single  time 
series  is  shown  in  Figure  5.1.  In  Figures  5.1  through  5.4,  the  normalized  frequency 
0.5  corresponds  to  10  KHz  because  the  sampling  rate,  /«,  was  20  KHz.  The  time 
series  were  multiplied  by  a  complex  exponential  of  the  form 
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Noraolizid  Tr^qumcy 

Figure  5.1:  Time  Series  Power  Spectra. 

so  as  to  shift  either  the  positive  or  the  negative  frequency  components  down  to  DC, 
depending  upon  the  sign  of  the  complex  exponential.  In  Equation  5.3,  /«  is  equal 
to  the  intermediate  frequency  of  5  KHz  as  discussed  in  Section  3.5.  The  spectrum 
of  the  basebanded  time  series  is  shown  in  Figure  5.2.  Only  the  energy  at  DC  is 
desired,  so  an  FIR  filter  was  used  to  remove  the  double  frequency  component.  The 
frequency  response  of  the  FIR  filter  is  shown  in  Figure  5.3  and  the  filtered  spectra 
is  shown  in  Figure  5.4.  The  envelope  is  obtained  by  t2tking  the  maignitude  of  the 
complex  time  series.  Since  the  low  pass  filter  removes  approximately  half  of  the 
power  by  filtering  out  the  double  frequency  term,  a  factor  of  3  dB  must  be  added. 
The  factor  to  convert  from  digital  units  to  volts  is  (409.6)'^  The  gain  shown  in 
Table  3.3  is  subtracted,  and  the  stave  sensitivity  applied  to  convert  from  voltage 
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to  pressure. 

5.4  Estimates  of  Bistatic  Surface  Scattering  Strengths 

Estimates  of  bistatic  surface  scattering  strengths  were  calculated  for  each 
run.  Only  time  periods  during  which  a  single,  well-defined  scattering  patch  was 
found  were  used.  These  time  periods  are  shown  in  Table  5.1. 

Scattering  strength  estimates  are  shown  in  Figures  5.5  through  5.13.  The 
incident  grazing  angle  is  approximately  constant  during  a  run,  varying  by  at  most 
3*.  Individual  panels  in  each  figure  pertain  to  a  different  scattered  grazing  angle. 
The  bistatic  angle  is  the  horizontal  axis  of  each  panel,  and  the  scattering  strength 
is  the  vertical  axis.  Table  5.2  shows  the  incident  scattering  angle,  t|>ine^  ^d  the 
range  of  scattered  grazing  angles  for  each  figure.  Note  that  for  fwd-las.108  and 
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Table  5.1:  Bistatic  Surface  Scattering  Strength  Time  Periods. 


Run 

Begin  Time  (msec) 

End  Time  (msec) 

fwd-las.102 

715 

800 

fwd-las.103 

724 

800 

fwd-las.104 

719 

800 

fwd-las.105 

718 

800 

fwd-las.108 

539 

550 

fwd-las.112  (30KHz) 

707 

750 

fwd-las.112  (40KHz) 

707 

750 

fwd-las.116 

532 

575 

fwd-las.117 

531 

575 

Table  5.2:  Range  of  Angles  Shown  in  Bistatic  Surface  Scattering  Strengths. 


Run 

t/’inc 

minV’.cof 

maxV'.cat 

fwd-las.102 

9.5 

10 

30 

fwd-las.103 

9.5 

10 

30 

fwd-las.104 

3.3 

20 

60 

fwd-las.105 

3.3 

20 

70 

fwd-las.108 

48.5 

6 

6 

fwd-las.112  (30KHz) 

9.3 

15 

35 

fwd-las.112  (40KHz) 

9.3 

15 

35 

fwd-las.116 

18 

11.5 

11.5 

fwd-las.117 

12.5 

20 

40 
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fwci-las.116  there  is  only  one  scattered  grazing  angle  and  therefore  only  one  panel 
in  these  figures.  While  the  display  of  this  information  in  panel  form  for  each  run 
may  appear  to  be  less  than  compact,  recall  that  bistatic  surface  scattering  strength 
depends  upon  three  geometrical  parameters,  ^»eat  &nd  as  well  as  wind 
speed  which  describes  the  surface  conditions.  Any  attempt  to  combine  the  results 
into  a  single,  three-dimensional  display  invariably  hides  information  which  must  be 
displayed.  In  Chapter  6  these  plots  are  again  shown  with  theoretical  predictions 
overlaid,  and  the  efficacy  of  this  display  format  will  be  evident. 

The  most  evident  trend  shown  in  the  estimates  of  the  bistatic  surface  scat¬ 
tering  strength  is  that,  over  the  narrow  range  of  bistatic  and  scattered  grazing 
angles  for  which  data  are  available,  the  scattering  strengths  are  fairly  constant. 
All  but  two  of  the  runs,  fwd-las.108  and  fwd-las.116,  have  a  bearing  angle  greater 
than  or  equal  to  55**,  which  causes  most  of  the  usable  time  periods  to  have  bistatic 
angles  less  than  —45*’.  Also,  the  wind  speed  is  at  least  4  m/s  for  most  of  the  runs, 
which  causes  much  of  the  usable  data  to  be  dominated  by  scattering  from  the  sub¬ 
surface  bubble  layer.  Derived  scattering  strengths  are  roughly  independent  of 
bistatic  angle  because  scatter  from  a  sub-surface  bubble  layer  is  roughly  isotropic. 

In  the  plots  of  runs  fwd-las.104,  105,  112-30KHz  and  117,  at  the  largest 
scattered  grazing  angle  and  the  smallest  bistatic  angles,  corresponding  to  the  region 
close  to  forward  scatter,  a  slight  increase  in  the  scattering  strength  can  be  seen. 
Here  the  effect  of  surface  roughness  is  becoming  dominant  over  the  effect  of  bubble 
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Scattered  Grazing  Pngle  =  6  deg. 


Bistatic  Angle,  deg. 


Figure  5.9:  Bistatic  Surface  Scattering  Strengths  fwd-las.108  (wind  speed  =  1.4 
m/s  and  V’ine  s  48.5”). 


scattering,  as  will  be  seen  in  Chapter  6  when  the  theoretical  predictions  are  overlaid 
on  the  plots  shown  in  this  section. 


5.5  Summary  and  Conclusions 


The  process  of  estimating  bistatic  surface  scattering  strengths  began  with 
a  test  for  homogeneity  to  avoid  processing  data  that  was  non-stationary.  With  the 
data  verified  as  stationary,  the  magnitude  of  the  complex  envelope  was  calculated 
using  the  method  described  in  Section  5.3.  The  calculation  of  the  bistatic  surface 
scattering  strength  was  done  using  the  sonar  equation  written  in  Equation  5.2. 
The  procedure  by  which  the  scattering  patch  area  and  location  was  determined  is 
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Figure  5.11:  Bistatic  Surface  Scattering  Strengths  fwd-las.112  (40  KHz)  (wind 
speed  3=  5.3  m/s  and  xlfi^e  =  9.3*). 
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Figure  5.12:  Bistatic  Surface  Scattering  Strengths  fwcl-las.116  (wind  speed  =  5.3 
m/s  and  =  18.0"). 

explained  in  Section  5.2,  and  the  computer  code  is  listed  in  Appendix  B.  While  the 
code  and  the  method  by  which  it  locates  the  scattering  patch  may  not  be  elegant, 
there  are  a  number  of  checks  included  to  avoid  choosing  an  erroneous  location  or 
size  for  the  scattering  patch.  The  horizontad  array  alone  only  is  used  in  order  to 
reduce  the  scattering  patch  area  ambiguity  as  detailed  in  Section  3.3.1. 

The  results  show  that  for  the  narrow  range  of  values  of  geometries  for  which 
data  are  available,  bistatic  scattering  strength  appears  to  be  roughly  independent 
of  all  three  geometrical  parameters.  This  results  from  most  of  the  data  being 
dominated  by  sub-surface  bubble  scatter,  which  itself  is  roughly  isotropic. 
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Ch  ipter  6 


Comparison  of  Bistatic  Surface  Scattering  Strength 
Estimates  with  the  McDaniel  Model 

6.1  Introduction 

Bistatic  surface  scattering  strength  estimates  from  the  FLIP  data  are  com¬ 
pared  with  bistatic  surface  scattering  strength  predictions  made  using  a  theory 
developed  S.T.  McDaniel  (15].  The  means  of  comparison  with  McDaniel’s  ex¬ 
pression  was  a  direct  comparison  of  the  bistatic  surface  scattering  strengths  shown 
in  Chapter  5  with  those  predicted  by  the  theory. 

6.2  S.T.  McDaniel’s  Expression 

The  expression  for  bistatic  surface  scattering  strength  developed  by  S.T. 
McDaniel  [15]  models  the  several  effects  that  may  be  present  in  scattering  from  the 
ocean  surface.  They  are  scattering  from  the  surface  itself,  which  is  modeled  by  the 
Kirchhoff  approximation  near  specular  and  by  the  Rayleigh- Rice  approximation 
far  from  specular;  isotropic  scattering  from  the  sub-surface  bubble  layer;  and  the 
attenuation  of  sound  traveling  to  and  from  the  surface  through  a  sub-surface  bubble 
layer.  Section  6.2  is  taken  wholly  from  [15]  and  is  discussed  here  in  detail  because 


[15]  may  not  be  easy  to  obtain. 

Scattering  from  the  surface  itself  is  divided  into  three  angular  regions: 

•  near  the  specular  direction  (KirchhofF  approximation), 

•  far  from  the  specular  direction  (Rayleigh-Rice  model), 

•  a  crossover  region  between  the  two  prior  regions  (weighted  sum  of  the  above). 

One  can  delineate  the  different  geometrical  regions  by  looking  at  the  normalized 
magnitude  of  the  sum  of  the  horizontal  components  of  the  incident  and  scattered 
wave  number  vectors.  It  is  bounded  by  the  relation 

0  <  <  cos  V».ne  +  COS  IpBcat  (6- 1 ) 

1^1 

where  itx  is  the  horizontal  component  of  the  incident  wave  number  vector  and  is 
the  horizontal  component  of  the  scattered  wave  number  vector.  The  wave  number 
vectors  of  the  incident  and  scattered  wave  are 

k  =  [-XCOSV'inc  +  2sini/;inc]fc  (6‘-) 

and 

k'  =  [f  cos  ipicat  cos  0  +  y  cos  V’»ca«  sin  0  +  zsin  ri’teatjk  (6.3) 

where  k  =u;/c.  The  horizontal  components  are  the  coefficients  of  x  and  y.  The  nor¬ 
malized  magnitude  of  the  sum  of  the  horizontal  components  can  then  be  expressed 
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in  terms  of  the  angles  as 


I  ,  2  2  1/1 

- I^r -  =  (cos*  1p,nc  +  cos*  lltscat  “  2  COS  4\nc  COS  V'jcol  COS  0)*'*. 

\k\ 


(6.4) 


In  the  specular  direction,  Equation  6.4  is  equal  to  0  because  and  0 

is  equal  to  0**.  In  the  backscattered  direction,  Equation  6.4  is  equal  to  2cosV’inc 
because  0,nc  =  i’tcat  and  0  is  equal  to  ISO**. 


6.2.1  Near-Specular  Scattering 

The  near  specular  region  can  be  defined  mathematically  \s 


0  <  (cos^  il’inc  +  COS^  “  2  COS  ti»me  COS  V’*co<  COS  .  (6.5) 


where  must  be  obtained  empirically.  In  this  region,  scattering  is  modeled  using 
the  KirchhofT  approximation  [13].  The  scattering  strength  can  then  be  expressed 
as 


S/{  —  lOlogjQ 


<  {n  kf  >  (-Q^  -I3^\ 

Alfflrfiy 


(6.6) 


where  a  is  the  x-component  of  the  sum  of  the  incident  and  scattered  wave  number 
vectors 


Q  =  A:(cOS  0,COJ  cos  <i>  —  cos  V’iTie)i 


(6.7) 
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0  is  the  y-component  of  the  sum  of  the  incident  and  scattered  wave  number  vectors 

0  =  kcos  i>tcat  sin  0,  (6.8) 

fix  is  the  product  of  the  z-component  of  the  sum  of  the  incident  and  scattered  wave 
number  vectors  and  the  x-component  of  the  large  scale  rms  surface  slope,  c*', 

fix  =  k{s\n  il>xcat  +  sin  (6.9) 

fly  is  the  product  of  the  z*component  of  the  sum  of  the  incident  and  scattered  wave 
number  vectors  and  the  y-component  of  the  large  scale  rms  surface  slope,  Cy, 

fly  =  fc(sin  +  sin  V»,nc)cry',  (6.10) 


and  n  is  the  outward  normal  to  the  surface.  The  normal  to  the  surface  can  be 
approximated  in  terms  of  the  local  surface  slopes  as 


n 


—z  — 


^  dx  dy 


(6.11) 


where  dCldx  and  dQjdy  are  the  local  surface  slopes  and  C(x,  y)  is  the  vertical  surface 
displacement.  A  more  physical  explanation  for  the  terms  in  the  exponential  in 
Equation  6.6  is  that  each  term  is  a  ratio  of  the  slope  at  the  point  (x,  y,  z)  ensonihed 
by  the  incident  wave  to  the  rms  surface  slope.  The  exponential  in  Equation  6.6 
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can  then  be  thought  of  as  the  probability  that  the  surface  will  be  correctly  tilted  to 
provide  a  specular  path  between  the  source  and  the  receiver  given  that  the  slopes 
are  Gaussian  distributed. 

At  sufficiently  high  frequency,  the  large  scale  surface  slope  is  independent 
of  direction  [15]  and  the  slope  components  can  be  expressed  in  terms  of  the  total 
large  scale  rms  slope  a\  i.e. 

=  (Ty' =  (6.12) 

Using  the  incident  wave  number  vector  shown  in  Equation  6.2  and  the  normal 
vector  approximation  shown  in  Equation  6.11,  the  expectation  in  Equation  6.6  can 
be  evaluated.  The  mean  slopes  are  0  and,  from  Equation  6.12,  the  second  moments 
are  a'^/2.  Mathematically,  this  is 

<  (^)  >=  0,  (6.13) 

<(|^)>=0,  (6.14) 

dy 

<(|^)’>=.7''/2,  (6.15) 

<(|)’>=.'’/2. 


(6.16) 


Therefore,  the  expectation  in  Equation  6.6  can  be  written  as 


<  (n  •  kf  >=  A:^(sin*  +  —  cos^  V'.nc) 


(6.17) 


Defining  the  angle  term  in  Equation  6.5  to  be  the  scattering  strength,  5^,  can 
be  expressed  as 


Sa  =  10  logio 


2  sin*  Ipinc  +  «■'*  cos*  xl^inc  1 

-0*  \\ 

4jr(sin  V'.nc  +  sin^,co«)*<T'^  \ 

,  (sin  y;,nc  +  sin  )*«^'^  / . 

(6.18) 


in  dB/rn^.  The  method  by  which  <t'^  is  found  is  described  in  Section  6.2.6. 


6.2.2  Far  From  Specular  Scattering 

Now  define  V'j  as  the  lower  bound  on  the  region  which  is  far  from  specular, 
which  can  be  modeled  using  the  Rayleigh-Rice  approximation  [12]. 


<  (cos^  +  cos*  ^,cot  -  2  cos  COS  r/»,cot  cos  (^)‘^*  <  2  (6.19) 


The  value  of  V'3  also  must  be  determined  empirically.  The  scattering  strength,  Sb, 
can  be  expressed  as 


Sb  =  10  log, 0 


<  (n-ky{n-P)^  > 

TT 


(6.20) 
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in  dB/m^  where  W  denotes  the  ocean  surface  wave  number  spectrum.  The  wave 
number  spectrum  is  omnidirectional  for  large  arguments,  so  that  using  Equa¬ 
tion  6.4,  it  can  be  written  as 


W{ki,  +  k^)  =  iy[A:{cOS^  v>.nc  +  COS^  -  2  cos  1p,nc  cos  0,cal  COS  0)*^*].  (6.21 ) 


The  expectation  in  Equation  6.20  can  be  evaluated  2is  in  Section  6.2.1  with  the 
additional  knowledge  that  for  a  Gaussian  slope  distribution 


<  >=  — 


(6.22) 


and  the  assumption  that  the  slope  components  are  uncorrelated 


(6.23) 


The  expectation  becomes 


<  (n  •  fc)’(n  •  k')^  >  =  k*  sin*  0,cat  + 

<T^(cos*0,cot  sin*  0inc  +  COS*  0>,nc  sin*  0,c,() 

2 


4<y^(3in  0>,nc  cos  0>inc  sin  0,cat  cos  cos  0) 
2 

<7**  cos*  t/>mc  COS*  0,cttt(3  COS*  0  sin*  0)' 

4 


-I- 


(6.24) 
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With  Equation  6.21  and  Equation  6.24,  values  for  Sb  can  be  calculated  using 
Equation  6.20.  The  method  for  calculating  the  surface  wave  number  spectrum  is 
discussed  in  Section  6.2.6. 

6.2.3  Crossover  Region  Scattering 

For  the  region  between  V’l  and  in  which  neither  the  KirchhofF  approxi¬ 
mation  nor  the  Rayleigh-Rice  approximation  can  be  exclusively  applied,  a  weighted 
sum  of  the  two  is  used  to  predict  the  scattering  strength.  This  crossover  region  is 
defined  mathematically  as 


V’l  <  (cos’  tl>inc  +  cos’  V'jc®*  —  2  cos  tine  COS  tseat  COS  <  ^2  (6.25) 

with  <02-  Defining  the  crossover  scattering  strength  as  5c,  the  weighted  sum 
is  written 

Sc  =  10 log,o(.4  max( +  .6(10^^''/*®’  -I- 10^^*'*®’))  (6.26) 

in  dB/rri^  where  Sa  and  Sb  are  defined  in  Sections  6.2.1  and  6.2.2.  These  three 
terms,  Sa,  Sb,  and  Sc  pertain  to  the  scattering  strength  of  the  surface  itself.  The 
effect  of  the  near-surface  bubble  layer  is  now  discussed. 
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6.2.4  Attenuation  Due  to  the  Bubble  Layer 

Two  effects  of  a  sub-surface  bubble  layer  are  modeled  in  S.T.  McDaniel’s 
expression,  those  being  attenuation  of  the  incident  and  scattered  rays  and  scattering 
of  sound  by  the  bubbles. 

Attenuation  of  the  incident  and  scattered  ray  paths  discussed  thus  far  by 
the  bubble  layer  is  modeled  as  lOlogfA),  where  A  is 


jmd 


A  =  exp 


+ 


Q 


sin  V»,ea» 


)l 


(6.27) 


2T'^ar^N 

6r 


(6.28) 


Ar,  expressed,  in  meters  is  the  resonant  bubble  radius  and  is  related  to  frequency 
by 

Or  =  3.2//  (6.29) 


where  /  is  the  transmit  frequency  in  Hz.  N  is  the  depth  integrated  bubble  density, 


N 


(6.30) 


where  N{z)  is  the  number  of  resonant  bubbles  per  unit  volume  at  a  depth  z. 
Expressions  for  as  a  function  of  wind  speed  and  the  resonant  bubble  radius  have 
been  developed  [14].  The  expression  for  N  is  also  different  depending  upon  whether 
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the  location  is  in  the  open  ocean  or  in  saline  bays/inlets.  Significantly  higher  bubble 
densities  have  been  found  in  saline  bays  and  inlets  [15).  The  reradiation  damping 
constant  is  dr  —  0.0136. 

The  expression  in  Equation  6.27  is  the  attenuation  at  a  particular  frequency, 
/,  and  integrated  bubble  density,  iV,  which  is  imparted  to  the  rays  which  travel  to 
and  from  the  surface.  As  either  the  incident  or  the  scattered  grazing  angle  goes 
to  0®,  implying  propagation  parallel  to  the  surface,  the  quantity  lOlog(A)  goes 
to  minus  infinity.  The  scattering  strength  calculated  in  Sections  6.2.1,  6.2.2  and 
6.2.3  are  reduced  by  a  factor  of  101og(i4)  to  account  for  attenuation  due  to  the 
sub-surface  bubble  layer. 

6.2.5  Scattering  Due  to  Bubble  Layer 

In  addition  to  the  attenuation  mechanism  due  to  bubbles  described  in  Sec¬ 
tion  6.2.4,  the  sub-surface  bubble  layer  also  scatters  sound  directly.  The  derivation 
of  the  scattering  strength  due  to  the  bubble  layer  is  simplified  by  the  fact  that 
bubbles  scatter  isotropically.  The  scattering  strength,  in  dB,  due  to  the  bubble 
layer  is  defined  to  be 

101og,o(A5')  (6.31) 

where  A  is  the  attenuation  factor  and  5^  is  the  scattering  cross  section  of  a  resonant 
bubble  of  radius  Or. 

The  attenuation  factor  is  the  sum  of  the  attenuations  for  the  four  possible 


82 


paths  from  the  source  to  the  bubble  scatterer  and  to  the  receiver: 

•  path  a:  direct  path  from  the  source  to  bubble  scatterer  to  the  receiver, 

•  path  b:  one  bounce  path  from  the  source  to  bubble  scatterer  to  the  surface 
and  to  the  receiver, 

•  path  c:  one  bounce  path  from  the  source  to  the  surface  to  bubble  scatterer 
and  to  the  receiver. 

•  path  d:  two  bounce  path  from  the  source  to  the  surface  to  bubble  scatterer  to 
the  surface  to  the  receiver. 


Note  that  the  four  paths  reduce  to  two  paths  when  the  geometry  is  monostatic.  To 
calculate  the  bistatic  scattering  strength  due  to  a  bubble  scatterer,  the  contribution 
and  attenuation  due  to  all  four  paths  must  be  included. 

The  attenuation  losses  are  expressed  in  the  same  form  as  in  Equation  6.27 
for  each  of  the  four  paths.  The  attenuation  for  the  four  paths  are  denoted  by 
where  n  and  m  indicate  the  path  [14].  The  attenuation  components  are: 

path  a: 

CiW  =  exp  [-  —  +  -r^ — \  F{2)  ,  (6.32) 

[  \sin^i„c  sinV'.cat/ 


path  b: 


2  2 
<^1^4 


f  (j) 

sin  V’.caJ  / 


2Q 

sin  tpscat 


(6.33) 
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path  c: 


2  2 

02  03  =  exp 


l-t 


^ - 

sin  Ptcat  sm  li’inc )  sin  Vine 


(6.34) 


path  d: 


V=exp[(j 


sin0i„c  smtf;tcatj  sin^,ca«  sm^inc 


(6.35) 


Here  F(s)  is  the  fraction  of  bubbles  below  the  depth  z,  which  can  be  expressed  in 


terms  of  N{2)  by 


dFjz)  _  N{z) 
dz  N 


(6.36) 


F(2)=^J"N(z')dz'. 


(6.37) 


Using  Equations  6.32  through  6.35,  the  attenuation  factor  is  equal  to 


^  =  [  N{z)an'^am‘^dz. 

^1^3  Jo 


(6.38) 


Equation  6.38  can  be  solved  for  each  term  in  the  summation  by  substituting  NdF{z) 
for  N{z)dz.  The  resulting  expression  is 


A  =  ^-[l-e 
<?  +  <?' 


-2(7+?')|  _  e-2»'j 

q'  -q 


(6.39) 
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where 


9  =  - 


sin 


(6.40) 


and 


9  =  - 


sin  Xl>tcat 


(6.41) 


and  Q  is  defined  in  Equation  6.28. 

The  scattering  cross  section  area  of  a  bubble  of  radius  Ur  in  the  layer  is 
equal  to 


52  = 


irOr 

~26 


(6.42) 


where  6  is  the  damping  factor,  which  is  related  to  the  reradiation  damping  constant 
Sr  by  the  relation 

6  =  6,+  .3.8  X  (6.43) 


Note  that  the  expression  shown  in  Equation  6.42  is  for  a  single  scatterer  of  radius 

а,  since  the  expression  for  the  attenuation  factor,  A,  was  integrated  over  the  entire 
water  column.  With  the  results  shown  in  Equations  6.39  and  6.42,  the  scattering 
strength  due  to  the  bubble  layer,  Sbub  c&n  be  calculated  using  Equations  6.31, 

б. 39,  and  6.42. 


6.2.6  Calculation  of  Large  Scale  Surface  Slope  (<7') 

The  quantities  <t'  and  W{K)  are  required  in  order  to  calculate  the  sur¬ 
face  scattering  strength.  The  two-dimensional  ocean  wave  height  spectra  W{I\)  is 
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related  to  the  one-dimension<il  wave  height  spectra  W{K)  via  the  relationship 

W{K)  =  (6.44) 

since  the  oce2ni  wave  height  spectrum  is  assumed  to  be  omnidirectional  at  high 
wavenumbers  as  was  discussed  in  Section  6.2.2.  References  [21],  [22],  and  [23] 
contain  appropriate  models  for  computing  W(I\)  as  a  function  of  wind  speed.  The 
large  scale  rms  surface  slope  c'  can  then  be  calculated  using 

K^W{K)dK  (6.45) 

Jo 

where  a  cutoff  value  of  A'l  =  0.4fc  has  been  found  to  be  satisfactory  [15]. 

6.3  Comparison  with  Bistatic  Surface  Scattering  Strengths 

Comparison  of  the  McDaniel  theory  predictions  with  the  results  presented 
in  Section  5.4  will  be  done  using  the  same  presentation  style.  The  results  shown  in 
Figures  5.5  through  5.13  will  now  be  shown  with  the  theory  predictions  overlaid. 
An  overview  of  the  data  presented  in  this  section  was  given  in  Table  5.1  and  5.2. 
The  data  were  described  in  more  detail  in  Section  5.4. 

One  preliminary  comment  is  made.  The  theory  prediction  is  largely  de¬ 
pendent  on  the  surface  wave  height  spectrum  to  describe  the  surface.  The  surface 
wave  height  spectrum  is  based  upon  wind  speed.  However,  the  surface  must  be  fully 
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Figure  6.1:  Wind  Speed  as  a  Function  of  Time. 

developed  before  the  wind  speed  provides  an  accurate  estimate  of  the  surface  wave 
height  spectrum.  To  be  fully  developed  requires  a  sustained  period  of  constant 
wind  speed  from  the  same  direction.  This  is  rarely  the  case  in  an  ocean  environ¬ 
ment,  so  that  the  appropriateness  of  the  calculated  surface  wave  height  spectrum 
is  always  an  open  question.  Figure  6.1  shows  the  wind  speed  as  a  function  of  time 
during  the  FLIP  experiment,  and  Figure  6.2,  the  wind  direction.  When  measured 
scattering  strengths  show  significant  deviation  from  the  theory  prediction,  it  can 
sometimes  be  explained  by  variable  wind  conditions  during  the  time  preceding  the 
measurement. 

The  predictions  of  the  theory  are  presented  overlaid  on  the  data  in  Fig¬ 
ures  6.3  through  6.11.  This  presentation  format  is  chosen  to  display  the  variation 
of  scattering  strength  with  bistatic  angle,  given  that  the  incident  amd  scattered 
grazing  angles  and  all  other  variables  are  fixed.  The  wind  speed,  frequency  and 
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Figure  6.2:  Wind  Direction  as  a  Function  of  Time. 

incident  grazing  angle  are  fixed  for  each  figure  or  run.  Each  panel  in  a  figure  per¬ 
tains  to  a  different  scattered  grazing  angle,  which  translates  to  a  different  time  after 
transmit  during  the  run.  Two  curves  for  the  prediction  are  shown.  The  line  at  a 
constant  value  is  the  scattering  strength  due  to  the  sub-surface  bubble  layer.  Its 
constant  value  with  bistatic  angle  is  due  to  the  isotropic  nature  of  scattering  from 
a  sub-surface  bubble  layer.  The  second  curve  is  a  combination  of  the  Kirchhoff 
approximation  and  the  Rayleigh- Rice  approximations  as  detailed  in  Sections  6.2.1 
through  6.2.3. 

Although  two  curves  are  presented  here,  in  ail  cases  the  theory's  prediction 
for  the  scattering  strength  would  be  the  higher  curve  at  a  particular  bistatic  angle. 
In  Figure  6.3  for  example,  at  a  scattered  grazing  angle  of  10*  (the  upper  left  panel) 
and  bistatic  angles  from  0*  to  -25*,  the  Kirchhoff/Rayleigh-Rice  curve  is  higher 
and  would  be  used.  At  bistatic  angles  from  —25*  to  —180*,  the  bubble  layer  curve 
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would  be  used.  The  two  curves  are  both  shown  for  all  bistatic  angles  to  illustrate 
where  each  mechanism  is  dominant. 

The  rate  of  decrease  in  the  scattering  strength  prediction  with  bistatic  angle 
near  the  =  0**  peak  increases  with  decreasing  wind  speed.  The  higher  the  wind 
speed,  the  lower  and  broader  the  peak  will  be  around  a  ^  =  0°,  which  corresponds 
to  forward  scattering. 

The  comparisons  between  runs  fwd-las.102  and  103  and  the  prediction, 
Figures  6.3  and  6.4,  show  rms  deviations  of  4.0  and  3.8  dB,  respectively.  The 
wind  speed  at  the  time  of  the  measurement,  7.2  m/s,  was  indicative  of  the  wind 
speeds  on  January  17‘^  and  the  agreement  with  the  theory  reflects  this.  All  of 
the  data  points  fall  in  the  region  which  the  theory  predicts  to  be  dominated  by 
bubble  layer  scattering.  A  few  of  the  panels  contains  no  data,  but  the  predictions 
are  included  for  comparison.  Runs  fwd>la5.104  and  105,  shown  in  Figures  6.5 
and  6.6,  exhibit  even  better  agreement  between  the  data  and  the  predictions,  with 
rms  deviations  of  2.0  and  1.1  dB  respectively.  The  wind  speeds  for  the  two  runs, 
5.5  m/s  and  5.8  m/s  respectively,  are  in  the  lower  range  of  wind  speeds  recorded 
on  the  day  the  measurements  were  made,  which  may  be  why  the  prediction  is  in 
general  lower  than  the  data.  The  data  in  both  runs,  for  the  most  part,  occur  in 
the  region  which  is  dominated  by  bubble  scattering.  A  few  points  at  high  scattered 
grazing  angles,  >  55®,  are  in  the  crossover  region  between  bubble  scattering  and 
Kirchhoflt/Rayleigh-Rice  scattering.  Figure  6.7  shows  the  comparison  for  run  fwd- 
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Figure  6.4:  Comparison  of  Bistatic  Surface  Scattering  Strengths  with  Theoretical 
Predictions  for  run  fwd-las.103  (wind  speed  =  7.2  m/s  and  r/li,nc  =  9.5"). 
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Figure  6.6:  Comparison  of  Bistatic  Surface  Scattering  Strengths  with  Theoretical 
Predictions  for  run  fwd-las.105  (wind  speed  =  5.8  m/s  and  ij^ine  =  3.3**). 
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Scattered  Grazing  fingle  =  6  deg. 


Figure  6.7:  Comparison  of  Bistatic  Surface  Scattering  Strengths  with  Theoretical 
Predictions  for  run  fwd-las.l08  (wind  speed  =  1.4  m/s  and  0,ne  *  48.5®). 

las.  108,  and  here  the  data  were  dominated  by  scattering  from  surface  itself.  This 
is  a  result  of  the  low  wind  speed,  1.4  m/s,  and  the  receiver  bearing  of  8*.  The  data 
match  the  theory  with  an  rms  deviation  of  5.4  dB.  The  wind  speeds  on  January 
21*‘  were  quite  low,  with  1.4  m/s  being  toward  the  middle  of  the  range  and  several 
measurements  as  low  a  1  m/s.  The  prediction  is  higher  than  the  data,  perhaps 
because  the  sea  surface  was  not  a  fully  developed  1.4  m/s  wind  driven  surface. 

The  multi'frequency  run,  fwd*las.ll2,  which  occurred  at  30  and  40  KHz, 
displayed  some  frequency  dependence  in  the  scattering  strengths.  The  theory  pre¬ 
dicts  only  a  slight  increase  in  bubble  layer  scattering  strength  as  the  frequency  is 
increased.  At  both  30  and  40  KHz  the  data  and  the  prediction  agreed  quite  well. 


Slrxgth.  (•  SoUtaring  Slr«9ll..  <*  i-oll.  .ng 
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5cctte'"ea  3rc2;>g  ^rgle  =  11.5  cec. 


Figure  6.10:  Comparison  of  Bistatic  Surface  Scattering  Strengths  with  Theoretical 
Predictions  for  run  fwd-las.116  KHz  (wind  speed  s  5.3  m/s  and  ~  18.0*). 

The  comparison  for  run  fwd*las.ll6  shows  good  agreement,  with  3.6  dB  rms 
deviation  between  the  data  and  the  theory.  Figure  6.10  is  similv  to  Figure  6.7,  with 
only  one  panel  showing  an  incident  grazing  angle  of  18*  and  a  scattered  grazing 
angle  of  11.5*.  The  data  are  clustered  in  a  region  which  the  theory  predicts  to  be 
dominated  by  bubble  layer  scattering,  and  the  relatively  constant  value  of  the  data 
for  small  changes  in  bistatic  angle  is  consistent  with  the  theory.  Run  fwd-las.117 
shows  good  agreement  between  the  data  and  the  theory.  Data  are  available  for  the 
regions  dominated  by  bubble  layer  scattering,  20.0  <  rl>tent  <  25.0,  and  by  surface 
roughness,  27.5  <  S  40.0.  Both  fwd-las.116  and  117  were  taken  within  15 
minutes  of  each  other  (see  Table  3.4  for  exact  times).  The  bearing  angle  was  30* 
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Table  6.1:  RMS  Deviations  between  the  Model  and  the  Data. 


Run 

RMS  Deviation  (dB) 

fwd-las.102 

4.0 

fwd-las.103 

3.8 

fwd-las.104 

2.0 

fwd-las.105 

1.1 

fwd-las.108 

5.4 

fwd-las.112  {30KHz) 

1.1 

fwd-las.112  (40KHz) 

1.6 

fwd-las.116 

3.6 

fwd-las.117 

2.9 

for  fwd*las.ll6  and  60"  for  fwd*las.ll7,  giving  a  different  set  of  grazing  angles  and 
bistatic  angles.  Table  6.1  summarizes  the  rms  deviations  between  the  theory  and 
the  data  for  all  the  runs. 

6.4  Conclusions 

Direct  comparison  of  bistatic  scattering  strength  estimates  derived  from 
FLIP  data  with  predictions  by  the  McDaniel  theory  showed  good  agreement  over  a 
range  of  geometries  and  wind  speeds.  The  theory’s  predictions  for  bubble  scattering 
appear  to  be  valid  over  the  ranges  shown  below: 

•  -105"  <  0  <  -35", 

•  3.3"  <  0.nc  <  18", 


•  6"  <  <  40", 


U«vin(SlrvBlh.  iC  Sealtar«n  Slr»i9ll.,  I*  SuHtoing  Si.mik, 
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Figure  6.11:  Comparison  of  Bis*;atic  Surface  Scattering  Strengths  with  Theoretical 
Predictions  for  run  fwd-las.117  KHz  (wind  speed  =  4.9  m/s  and  0,„c  =  12.5*). 
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•  1.4  m/s  <  wind  speed  <  7.2  m/s, 

•  /  =  30,40  KHz. 

This  is  an  important  result  because  it  represents  the  first  such  validation  of  the 
McDaniel  model  prediction. 

Some  of  the  poorer  comparisons  may  be  due  to  an  improper  wind  speed 
used  to  characterize  the  surface  conditions.  The  theory  predicts  a  slight  increase  in 
bubble  scattering  strength  with  frequency  over  the  30  KHz  to  40  KHz  range.  The 
30  KHz  and  40  KHz  data  display  this  very  slight  trend  and  match  the  theory  well. 


Chapter  7 


Summary  and  Conclusions 


7.1  Introduction 

The  major  goal  of  the  research  presented  in  this  thesis  was  to  analyze  the 
statistical  characteristics  of  bistatic  surface  scattering  data  collected  during  the 
FLIP  experiment  of  January  1992,  and  to  compare  bistatic  scattering  strength 
estimates  to  the  McDaniel  prediction  model.  The  original  contributions  presented 
in  this  thesis  are  as  follows. 

•  The  ocean  surface  scatter  data  were  histogrammed  and  their  cumulative  distri¬ 
bution  function  compared  with  a  Rayleigh  distribution  using  the  probability  of 
false  alarm,  PpA-  Theory  predicts  that  the  dominant  surface  scattering  mech¬ 
anism  is  the  slope  of  surface  for  higher  grazing  angles  and  isotropic  resonant 
bubble  scatter  at  lower  grazing  angles  and  higher  wind  speeds.  The  magni¬ 
tude  of  the  complex  envelope  of  the  ocean  surface  scatter  data  was  found  to 
be  Rayleigh  for  both  mechanisms  and  for  all  but  the  highest  wind  speed  (7.2 
m/s).  This  result  shows  that  the  use  of  the  central  limit  theorem  to  describe 
the  statistical  characteristics  of  bistatic  surface  scatter  is  valid  over  the  geome¬ 
tries  and  wind  speeds,  except  7.2  m/s,  studied  here.  It  also  shows  that  the 


central  limit  theorem  is  valid  for  both  scattering  mechanisms,  surfa^re  rough¬ 
ness  and  the  bubble  layer,  for  the  geometries  and  all  but  the  highest  wind 
speed  studied  here. 

•  Mean  bistatic  surface  scattering  strength  was  calculated  using  the  sonar  equa¬ 
tion  and  compared  to  a  prediction  by  S.T,  McDaniel’s  theory  for  bistatic  sur¬ 
face  scattering  strengths  (15).  The  portion  of  the  theory  for  predicting  scatter¬ 
ing  strengths  dominated  by  bubble  scattering  was  validated  by  the  data  for  a 
number  of  geometries.  This  is  the  first  validation  of  that  portion  of  the  theory. 

•  A  technique  was  developed  for  combining  ocean  surface  scatter  data  for  differ¬ 
ent  grazing  angles  which  increased  the  number  of  data  points  in  the  ensemble. 
With  this  technique,  data  which  were  the  result  of  the  same  dominant  scatter¬ 
ing  mechanism,  the  surface  itself  or  the  sub-surface  bubble  layer,  were  com¬ 
bined  to  estimate  the  statistical  characteristics  of  the  scattering  mechanism. 

A  more  detailed  description  of  these  contributions  follows. 

7.2  Statistical  Characteristics  of  the  Data 

The  statistical  testing  was  comprised  of  tests  for  homogeneity,  normality 
and  for  Rayleigh  behavior.  Homogeneity  was  investigated  using  the  Kolmogorov- 
Smirnov  two  sample  test,  in  which  the  statistics  of  the  first  half  of  an  ensemble 
are  compared  with  those  of  the  second  half.  All  usable  data  were  found  to  be  sta¬ 
tistically  stationary.  Testing  for  homogeneity  is  a  necessary  first  step  in  the  data 
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analysis  since  a  change  in  the  underlying  statistical  characteristics  during  a  par¬ 
ticular  run  (or  ensemble)  would  indicate  a  meaningful  change  in  the  environment 
during  the  run.  In  this  case,  one  can  no  longer  assume  that  the  statistical  charac¬ 
teristics  of  the  entire  ensemble  come  from  the  same  parent  distribution,  rendering 
the  ensemble  invalid.  If  the  data  are  validated  as  homogeneous,  then  one  cam  treat 
each  return  in  the  run  as  an  independent  realization  of  the  reverberant  process. 

The  surface  scatter  time  series  were  tested  for  normality  using  the  two 
sample  Kolmogorov- Smirnov  test,  in  which  the  cumulative  distribution  of  the  data 
is  compared  to  that  of  a  Gaussian.  All  of  the  data  passed  the  K-S  test  for  normality. 

■  However,  deviations  from  normality  in  surface  scatter  reverberation  are  of  most 

I 

I 

interest  in  the  tails  of  the  distribution,  and  the  limited  number  of  the  pulses  (25 
to  50)  made  such  deviations  unobservable.  Also,  the  usefulness  of  the  K-S  test  to 
identify  non-normal  behavior  in  the  tails  of  the  distribution  is  questionable. 

In  order  to  increase  the  number  of  samples  used  in  the  test  for  normality, 
data  which  were  a  result  of  the  same  dominant  scattering  mechanism  were  nor¬ 
malized  and  grouped  together  into  a  single  ensemble.  This  combinational  method, 
which  is  described  in  Section  4.3,  increased  the  amount  of  data  in  all  runs  by  at 
least  a  factor  of  five,  resulting  in  better  definition  of  the  distribution.  The  longer 
ensembles  were  tested  for  normality  by  comparing  the  statistics  of  their  magnitude 
to  a  Rayleigh  distribution,  since  theory  predicts  that  the  magnitude  of  a  complex 
Gaussian  time  series  will  be  Rayleigh  distributed.  The  comparison  was  made  using 
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the  probability  of  false  alarm,  Pfa^  which  is  one  minus  the  cumulative  distribution, 
in  order  to  focus  on  the  largest  excursions  of  the  time  series  which  define  the  shape 
of  the  distribution  tail.  The  Pfa  comparison  showed  some  deviation  from  Rayleigh 
behavior  for  all  runs.  However,  only  at  the  highest  wind  speeds  of  7.2  m/s  were  the 
deviations  from  Rayleigh  behavior  significant.  While  apparent  1  ufficient  to 

fail  the  Kolmogorov-Smirnov  test,  the  deviations  do  indicate  that  at  higher  wind 
speeds,  the  distribution  tails  are  higher  than  those  of  a  Gaussian  distribution. 

7.3  Bistatic  Surface  Scattering  Strengths 

The  FLIP  bistatic  data  were  processed  to  yield  estimates  of  bistatic  surface 
scattering  strength  using  a  rather  inelegant  method  of  processing  as  detailed  in 
Section  5.2.  An  isovelocity  sound  speed  profile  was  used  because  ray  tracing  was 
unnecessary.  The  scattering  strength  estimates  are  shown  in  Figures  5.5  through 
5.13  as  a  function  of  scattered  grazing  angle,  a^nd  bistatic  angle  0,  for  each 
run,  which  in  turn  determines  wind  speed  and  incident  grazing  angle,  ^,„c.  No 
satisfactory  means  exists  to  combine  data  for  different  geometrical  p2U’ameters. 
The  majority  of  the  results  are  dominated  by  bubble  scattering,  as  can  be  seen  by 
the  their  very  weak  dependence  on  bistatic  angle.  This  characteristic  is  due  to  the 
isotropic  nature  of  scatter  from  bubbles. 
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7.4  Comparison  of  Scattering  Strength  Estimates  with  the  McDaniel 
Prediction  Model 

Bistatic  surface  scattering  strength  estimates  were  compared  with  predica¬ 
tions  of  a  theory  developed  by  S.T.  McDaniel,  which  is  presented  in  Section  6.2.  The 
theory  consists  of  two  main  parts,  one  modeling  scattering  from  the  surface  itself 
and  one  modeling  scattering  from  a  sub-surface  bubble  layer.  Scattering  from  the 
surface  is  modeled  by  either  one  or  a  combination  of  the  Kirchhoff  approximation 
and  the  Rayleigh-Rice  approximation.  Attenuation  by  the  sub-surface  bubble  layer 
of  scattering  from  the  surface  is  included.  Scattering  from  a  sub-surface  bubble 
layer  is  modeled  as  the  product  of  the  scattering  cross  section  of  a  bubble  resonant 
at  the  incident  frequency  and  the  depth  integrated  attenuation  factor.  The  depth 
integrated  attenuation  factor  accounts  for  four  paths  from  the  source  to  the  scat¬ 
tering  bubble  and  on  to  the  receiver,  and  for  the  bubble  density  integrated  over  the 
entire  water  column. 

The  comparisons  between  the  data  and  the  theory  showed  good  agreement, 
with  a  rms  deviation  of  about  3  dB  for  all  the  data.  Some  of  the  poorer  agreements 
may  be  a  result  of  the  wind  speed  being  much  higher  or  lower  during  the  period 
prior  to  the  run,  inasmuch  as  the  theory  uses  wind  speed  to  calculate  a  wave  height 
spectrum  which  has  a  strong  effect  on  the  prediction.  In  order  for  this  to  be  an 
accurate  and  valid  method,  the  wind  speed  must  have  been  fairly  constant  in  speed 


and  direction  for  some  time. 
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7.5  Recommendations  for  Future  Research 

Any  future  experiment  to  measure  bistatic  surface  scattering  would  benefit 
by  using  a  horizontal  array  only  to  receive  data  and  recording  each  hydrophone 
rather  than  only  the  staves.  This  would  increase  the  range  of  angles, 
and  for  which  statistics  could  be  generated,  since  the  beam  could  be  steered  over 
a  range  of  bearing  angles  instead  of  just  looking  perpendicular  to  the  array  axis  as 
was  done  with  the  data  collected  in  January  of  1992. 

From  an  analysis  point  of  view,  the  ability  to  steer  the  receive  beam  would 
permit  a  single  run  to  provide  estimates  of  the  bistatic  surface  scattering  strength 
over  a  large  range  of  bearing  angles  instead  of  just  a  few.  Two  or  three  runs  at 
bearing  angles  of  15®,  45®  and  75®  could  provide  data  for  all  bearing  angles  from  0® 
to  90®.  This  would  permit  the  study  of  some  of  the  more  subtle  angular  dependences 
of  bistatic  surface  scattering  strength.  While  more  computationally  intensive,  this 
would  permit  a  more  detailed  analysis  of  bistatic  surface  scattering  strengths  as  a 
function  of  geometry. 

Element  level  data  would  also  enhance  the  statistical  study  of  the  bistatic 
surface  scatter  since  beam  steering  could  be  used  to  look  at  the  surface  in  direc¬ 
tions  that  may  generate  non-Gaussian  scatter.  Perhaps  a  clearer  understanding  of 
the  transition  of  bistatic  surface  scatter  from  Gaussian  behavior  to  non-Gaussian 
behavior  would  emerge.  Increasing  the  number  of  pulses  transmitted  in  a  run 
would  also  increase  the  accuracy  of  the  statistical  analysis,  but  one  must  be  wary 
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of  non-stationarity  when  the  duration  of  the  runs  is  increased. 
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Appendix  A 

Fundamentals  of  Monostatic  and  Bistatic  Active  Sonar 

Boundary  Scatter 

In  order  to  understand  the  geometry  of  a  bistatic  active  sonar  boundary 
scatter,  one  should  begin  with  the  simpler  case  of  a  monostatic  active  sonar.  As¬ 
suming  an  isovelocity,  range- independent  sound  profile,  propagation  occurs  in  a 
straight  line. 

For  a  monostatic  sonar,  where  the  source  and  receiver  are  collocated,  and 
with  a  relatively  short  transmit  pulse  and  no  boundary  interaction,  the  ensonified 
volume  at  any  particular  time  is  a  spherical  shell,  the  radius  of  which  is  determined 
by  the  elapsed  time  since  pulse  transmission,  and  the  thickness  of  which  is  deter¬ 
mined  by  the  transmit  pulse  duration.  This  implies  that  at  a  particular  listen  time, 
only  volume  inhomogeneities  and  scatterers  located  in  this  shell  will  contribute  to 
the  return  seen  at  the  receiver. 

In  the  case  where  the  pulse  interacts  with  a  boundary,  such  as  the  ocean 
bottom  or  surface,  the  spherical  shell  is  cut  by  a  plane,  yielding  an  annular  ring 
on  the  boundary  which  contributes  to  the  return  at  particular  listen  time.  If  the 
received  signal  is  dominated  by  the  interaction  with  the  boundary,  one  can  ignore 
the  contribution  from  the  volume.  An  example  of  boundary  interaction  is  shown  in 
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Figure  A.l;  Overhead  View  of  Monostatic  Geometry. 

Figure  A.l.  We  further  define  the  scattering  patch  to  be  the  part  of  the  ensonified 
area  whose  contribution  to  the  return  is  within  3  dB  of  the  maximum  return  pro¬ 
duced  by  any  point.  The  directionality  of  the  source  and  receiver  determine  what 
portion(s)  of  the  ensonified  area  compromises  the  scattering  patch.  When  both 
the  source  and  the  receiver  are  omnidirectional,  the  scattering  patch  is  composed 
of  the  entire  ensonihed  area.  If  the  source  or  the  receiver  or  both  have  some  di¬ 
rectionality,  the  scattering  patch  will  be  some  smaller  portion(s)  of  the  ensonihed 
area.  Such  a  scattering  patch  is  shown  in  Figure  A.l  by  the  small  square  on  the 
x-axis.  Figure  A. 2  shows  a  side  view  of  the  monostatic  geometry  in  order  to  define 


Figure  A. 2:  Side  View  of  Monostatic  Geometry. 


the  incident  grazing  angle, 

In  the  case  of  a  bistatic  sonar  geometry,  the  source  and  the  receiver  are 
iorated  at  two  different  locations.  When  a  pulse  is  transmitted  by  the  source,  a 
delay  equal  to  the  length  of  time  for  the  pulse  to  propagate  from  the  source  to  the 
receiver  occurs  before  the  pulse  is  received  at  the  receiver.  Arrival  of  the  pulse  at  the 
receiver  via  a  direct  path,  without  any  boundau^y  interactions,  is  called  the  direct 
blast  in  the  vernacular  of  bistatic  geometries  [11].  After  the  direct  blast  ends,  only 
returns  from  scattering  due  to  inhomogeneities,  boundaries  and  scatterers  located 
in  the  ensonihed  area  are  received.  The  ensonihed  area  is  no  longer  a  spherical 
shell,  but  is  now  an  ellipsoidal  shell. 

When  the  ensonifted  volume  interacts  with  a  boundary,  the  resulting  en- 
sonified  area  is  an  elliptical  annulus.  An  example  of  the  ensonihed  area  on  a  surface 
is  shown  in  Figure  A. 3  looking  down  on  the  surface.  Three  geometrical  parameters 
are  needed  to  a  describe  bistatic  geometry,  as  shown  in  Figures  A.4  and  A. 5.  The 
bistatic  angle,  is  defined  to  be  zero  in  the  forward  scattered  direction.  The 

incident  grazing  angle,  V'inet  and  the  scattered  grazing  angle,  are  in  general 
not  equal,  as  shown  in  Figure  A.5. 


Figure  A.3:  Overhead  View  of  Bistatic  Geometry. 


Figure  A.4:  Definition  of  Bistatic  Angle  from  Top  View. 
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Appendix  B 

Code  to  Reduce  Data  to  Scattering  Strength 

B.l  Usage  of  Code 

The  code  included  in  this  appendix  was  used  to  reduce  the  bistatic  sur¬ 
face  scatter  data  collected  from  the  FLIP  experiment  an  produce  estimates  of  the 
bistatic  surface  scattering  strength. 

The  code  requires  the  user  to  input  a  number  of  values  unique  to  the  ge¬ 
ometry  and  operating  parameters  at  which  the  data  were  collected.  In  addition, 
input  files  are  required  for  the  mean  magnitude  of  the  complex  envelope  and  any 
directional  beam  patterns  used.  The  user  is  prompted  for  input  from  the  keyboard. 
The  output  consists  of  two  files,  both  in  ascii  format,  the  names  of  which  are  chosen 
by  the  user.  One  file  contains  the  estimate  of  bistatic  surface  scattering  strength 
as  a  function  of  the  three  geometrical  parameters,  ^,„e,  0«ea(  and  time.  The 

second  file  contains  the  position  of  the  center  of  the  scattering  patch  as  a  function 
of  time. 

B.2  Listing  of  Code 

This  is  a  listing  of  the  program  used  to  convert  the  data  from  an  estimate  of 
the  received  level  in  dB  into  an  estimate  of  the  bistatic  surface  scattering  strength. 
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37  C 

38  C 

39  C 

30  C 

31  C 
33  C 

33  C 

34  C 

35  C 

36  C 

37  C 

38  C 

39  C 

40  C 

41  C 
43  C 

43  C 

44  C 

45  C 

46  C 

47  C 

48  C 

49  C 

50  C 

51  C 
53  C 

53  C 

54  C 

55  C 

56  C 

57  C 

58  C 

59  C 

60  C 

61  C 
63  C 

63  C 

64  C 

65  C 

66  C 

67  C 

68  C 

69  C 

70  C 


progrw  v«r83 

Prograa  to  coaputo  th«  biatatic  aurfaca  acattaring  atrangth  for  tha 
FLIP  iata.  Thia  varaion  alloaa  for  both  tranaaittar  and  racaivar 
baaa  pattama  (tabla  drivan)  aith  arbitrary  point  diraction,  axcapt 
for  roll. 


Variablaa  uaad  in  prograa: 
Environaantal  Inforsation: 


c 

attn 


aound  apaad  (iaovalocity  in  a/a) 
attanuation  coafficiant  (function  of  fraq.) 
in  dB/a. 


Gaoaatry  and  Pulaa 


Information: 


trlaval 

duration 

f 

aaparation 


tranaaittad  laval  (dB) 
pulaa  langth  antarad  in  aaac  convartad 
to  aac 

valua  to  contort  bataaan  aaconda  and  aaac 
horizontal  diatanca  bataaan  tha  aourca 
and  tha  racaivar  (a) 


Tranaait  Baaa  Information: 


tr.baaa 

fila.tr 

baan_tr(32761) 

oriantation.tr 


alaa_fix.tr 


flag  for  tranaait  baaa  (  1  >  oani,  2  ■  tabla) 
fila  containing  tha  tranaait  baaa  pattam 
tranaait  baaa  pattam  atorad  in  array 
angla  of  HRA  for  tha  tranaaittar  aith  raapact 
to  lino  draan  bataaan  aourca  and 
racaivar  (CW  ia  *) 
angla  of  alavation  of  MU  for 
tranaaittar  (up  ia  *) 


Racoiva  Baaa  Information: 


rac.baam 

fila.rac 
baan_rac (32761) 
oriantation.rac 


alav.fiz.rac 


for  racoiva  baaa  (  1  ■  oani.  2  •  tabla) 
fila  containing  tha  racoiva  baaa  pattam 
racoiva  baaa  pattam  atorad  in  array 
angla  of  MU  for  tha  racaivarar  with  raapact 
to  lino  draan  bataaan  aourca  and 
racaivar  (CV  ia  '*') 

angla  of  alavation  of  MU  for  racaivar  (up  ia  *■) 


Envalopa  Infomation: 

fila.anv  fila  containing  tha  anvalopa 

raclaval(lOOO)  anvalopa  laval  of  raal  data  aamplod  at  Imaac 


Output  Fila  Infomation: 

fila.ottt  fila  containing  tha  paak  pointa  for 

aach  tima  alica 

fila.aa  fila  containing  tha  acattaring  atrangth  and  tha 

anglaa  for  aach  tima  alice 

Saarch  Infomation: 


atart 

and 

intarval 

naah.original 

maah 


maah.y 


atart  tima  for  calculation  of  allipaa 
and  tima  for  calculation  of  allipaa 
intarval  at  which  calculationa  ara 
to  ba  dona 

original  maah  aiza  input  by  uaar 
aiza  of  atop  takan  in  aaarch  for  pointa 
bataaan  tha  too  allipaaa  on  tha  aurfaca 
in  tha  x-diraction 

aiza  of  atop  takan  in  aaarch  for  pointa 
bataaan  tha  too  allipaaa  on  tha  aurfaca 
in  tha  y-diraction 
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ri  c 

72  c 

73  C 

74  C 

75  C 

76  C 

77  C 

78  C 

79  C 

80  C 

81  C 

82  C 

83  C 

84  C 

85  C 

86  C 

87  C 

88  C 

89  C 

90  C 

91  C 

92  C 

93  C 

94  C 

95  C 

96  C 

97  C 

98  C 

99  C 

100  C 

101  C 

102  C 

103  C 

104  C 

105  C 

106  C 

107  C 

108  C 

109  C 

no  c 

111  c 

112  C 

113  C 

114  C 

115  C 

116  C 

117  C 

118  C 

119  C 

120  C 

121  C 

122  C 

123  C 

124  C 

125  C 

126  C 

127  C 

128  C 

129  C 

130  C 

131  C 

132  C 

133  C 

134  C 

135  C 

136  C 

137  C 

138  C 


(liat.2_aurl 

(liat.2_rac 

diat 

ti«a 

round.trip.outar: 
rotmd.tr  ip.  innar : 


diatanca  from  aourca  to  aurfaca  point 
diatanca  Iroa  aurfaca  point  to  racaivar 
aua  of  diat_2.aurf  and  diat_2_roc 
ti»a  at  ahich  tha  calculation  ia  baing  dona 
total  traval  diatanca  to  tha  baginning 
of  tha  pulaa 

total  traval  diatanca  to  tha  and  of  tha  pulaa 


Coordinataa  in  un-tranalatad,  un-rotatad  cartaaian  coordinataa 


aourca.x  x  location  of  aourca 

aourca.y  y  location  of  aourca 

aourca.x  z  location  of  aourca  (down  ia  *) 

racaivar.x  x  location  of  racaivar 

racaivar.y  y  location  of  racaivar 

racaivar.x  x  location  of  racaivar  (dom  ia  '•’) 

dalta.a( 100000)  :  horizontal  diatanca  fron  aourca  to  racaivar 

dalta.bC 100000)  :  horizontal  diatanca  fron  aourca  to  allipaa  pt. 

dalta.c( 100000)  :  horizontal  diatanca  fron  racaivar  to  allipaa  pt. 

allipaa.xC 100000):  allipaa  x-coord. 

allipaa. y( 100000):  allipaa  y-coord. 

Coordiantaa  tranalatad  and  rotatad  for  aourca  oriantation 

aourca.x.tr  x  location  of  aourca  (aaparation) 

aourca.y.tr  y  location  of  aourca  (0.0) 

racaivar.x.tr  x  location  of  racaivar  (0.0) 

racaivar.y.tr  y  location  of  racaivar  (0.0) 

dalta.a( 100000)  :  horizontal  diatanca  froa  aourca  to  racaivar 

dalta.b( 100000)  :  horizontal  diatanca  froa  aourca  to  allipaa  pt. 

dalta.c ( 100000)  :  horizontal  diatanca  froa  racaivar  to  allipaa  pt. 

bmg.tr (100000)  :  baaring  i^la  froa  aourca  to  allipaa  pt. 

ralativa  to  HRA  of  aourca 

alav.tr (100000)  :  alavation  angla  froa  aourca  to  allipaa  pt. 

ralativa  to  HIU  of  aourca 

allip8a_x.tr( 100000) :  allipaa  z-coord.  rafarancad  to  aourca 
allipaa.y_tr( 100000) :  allipaa  y-coord.  rafarancad  to  aourca 
allipaa.z.tr( 100000) :  allipaa  z-coord.  rafarancad  to  aourca 

Coordiantaa  tranalatad  and  rotatad  for  racaivar  oriantation 


aourca.x.rac 
aourca.y.rac 
racaivar.x.rac 
racaivar.y.rac 
dalta.a( 100000)  : 
dalta.b( 100000)  : 
dalta.c (100000)  ; 
bmg.rac(  100000) : 

alav.rac ( 100000 ) : 

allipaa.x.rac( 100000) : 
allipaa.y.racdOOOOO) : 
allipaa.z.rac( 100000) : 


X  location  of  aourca  (0.0) 
y  location  of  aourca  (0.0) 

X  location  of  racaivar  (-laaaparation) 
y  location  of  racaivar  (0.0) 
horizontal  diatanca  fron  aourca  to  racaivar 
horizontal  diatanca  fron  aourca  to  allipaa  pt. 
horizontal  diatanca  fron  racaivar  to  allipaa  pt. 
baaring  angla  fron  racaivar  to  allipaa  pt. 

ralativa  to  lUU  of  racaivar 
alavation  angla  froa  racaivar  to  allipaa  pt. 

ralativa  to  HRA  of  racaivar 
allipaa  x-coord.  rafarancad  to  racaivar 
allipaa  y-coord.  rafarancad  to  racaivar 
allipaa  z-coord.  rafarancad  to  racaivar 


Saarch  for  Haxinun  Infomation: 


pointer 

pointar.rac 

nax 


location 
taat( 100000) 


pointer  for  location  in  trananit  bean  pattern 
pointer  for  location  in  racaiva  bean  pattern 
teat  variable  uaad  to  find  nazinun  in 
the  aun  of  the  trananiaaion  loaaaa. 
bean  pattema  and  attenuation  loaa 
location  of  tha  point  nax  in  tha  array  of  pointa 
batvaan  tha  tvo  allipaaa 
valuaa  corraaponding  to  tha  aun  of  tha 
trananiaaion  loaaaa.  bean  pattama  and 
tha  attenuation  loaa  for  each  point 
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139  C 

140  C 

141  C 
143  C 

143  C 

144  C 

145  C 

146  C 

147  C 

148  C 

149  C 

150  C 

151  C 

152  C 

153  C 

154  C 

155  C 

156  C 

157  C 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177  C 

178  C 

179  C 

180 

181  C 

182 

183  C 

184 

185  C 

186 

187  C 

188  C 

189  C 

190 

191 
193  C 

193 

194 

195 

196  C 

197 

198 

199  C 
300 
201 
303  C 

303 

304 

305  C 
206 
207 


Calculation  of  Aroa  Inforaation: 


aroa 


araa.f ino 

araa.finor 

flag 


aroa  of  onsonofiod  patch  in  a**2  onco  tho 
aroa  convorgos  to  a  valuos  (♦/-  lOX  of 

aroa  with  coaraor  nosh  aizo) 

aroa  oith  coarsor  nosh  aizo 
aroa  aith  anallor  noah  aizo 
flag  uaod  in  looping  to  finor  noah  aizo 


Roaulta  Information: 


biat.anglo 

thota_naught 

thota 

acat.atr 


biatatic  angla  (dog) 

incidont  anglo  from  aourco  to  onaonofiod  patch 
on  tho  aurfaco 

acattorod  anglo  from  onaonofiod  patch  on  tho 
aurfaco  to  tho  rocoivor 
valuo  of  acattoring  atrongth 


implicit  roaloS  (a-hi  k.  o-z) 

dimonaion  roclovol(lOOO) ,  ollipao_z( 100000) ,  ollipao.yC 100000) 
dimonaion  ollipao.z.rocdOOOOO) .  ollipao_y_roc(  100000) 
dimonaion  ollipao.z.rocdOOOOO),  ollipao.z.tr (100000) 
dimonaion  ollipao.y.trdOOOOO) ,  ollipao.z.trd00000) 
dimonaion  bmg.tr (100000 ) ,  olov.tr (100000) 
dimonaion  brag.roc (100000) .  olov.rac (100000) 
dimonaion  boam.tr (32761).  boam.roc(32761) 
dimonaion  dolta.a( lOOOOOK  dolta.bdOOOOO) ,  dolta.c (100000) 
dimonaion  toat (100000) 

roal*8  max.  moah.  moah.y.  x,  y.  moah.original 
intogor  countor,  atart.  and,  intorval,  flag 
roal  oriontation.tr.  olov.fix.tr 
roal  oriontation.roct  olov.fix.roc 
intogor  tr.boam,  roc.boam 
intogor  pointor.  pointor.roc 
intogor  poaition.tr.  poaition.roc 
intogor  bogin.  longth 

charactor«20  filo.tr.  filo.roc.  filo.onv,  filo.out.  filo.aa 


Pro~dofinod  valuoa  uaod  in  program: 

Sound  Spood:  c  (m/a) 

c  ■  1500.0 

Samplo  Rato:  f  (Hz) 

f  -  1000.0 

Convoraion  from  radians  to  dogrooa:  rad  (dog/rad) 

rad  •  57.2958 
Pi.  of  courao;  pi 

pi  -  3.141592654 

Asks  for  inputs  for  tho  run. 

«rito(o.o)  ’Entor  tho  transmittod  lovol  in  dB. ’ 
rosd(o,o)  trlovol 

orito(o.«)  'Entor  tho  pulso  duration  in  msoc. ' 
road(o,o)  duration 
duration  ■  durat ion/ 1000. 0 

vrito(o.o)  'Entor  tho  attouation  coofficiont  (dB/m).’ 
road(«.o)  attn 

vrito(o.o)  'Entor  tho  bogin  tiso  for  onvolopo  data  in  msoc. ' 
road(o,«)  bogin 

srito(o.«)  'Entor  tho  samplo  timo  for  onvolopo  data  in  msoc.’ 
road(o.o)  longth 

arito(o.o)  'Entor  tho  start  tins,  ond  tins,  and  intorval 
1  for  tho  run  in  msoc . ' 


US 


208 
209  C 
2t0 
211 
212 

213  C 

214 

215 

216  C 

217  C 

218  C 

219 

220 
221 
222 

223  C 

224 

225 

226  C 

227  C 

228  400 

229  1 

230  2 

231  C 

232  C 

233  C 

234  C 

235 

236  C 

237 

238 

239  1 

240 

241  C 

242 

243  1 

244 

245  C 

246 

247  C 

248 

249 

250  C 

251 

252  C 

253 

254  C 

255 

256  500 

257  1 

258  2 

259 

260  C 

261 
262 

263  1 

264 

265  C 

266 

267  1 

268 

269  C 

270 

271  C 

272 

273 

274  C 

275 

276  C 


r«ad(«,«)  start,  and,  intarvsl 

•fitsC*,*)  'Entsr  th«  assh  sizs  for  tho  surfscs  sosrch  in  nstors.  ’ 
rsad(«,«)  nosh 
Msh.originsl  >  nsah 

writoC*,*)  ’Entor  ths  sourcs-rscoivsr  sspsration. ’ 
rssd(«.*)  sopsrstion 

Pisco  tho  sourco  and  rocoivor  on  tho  z-axis  sywotrical  oith  oach  othor 

sourco.s  ■  -l«aoparat ion/2.0 
rocoivor _K  ■  soparat ion/2.0 
sourco.y  >0.0 
rocoivor.y  >0.0 

srito(o,«)  ‘Entor  tho  sourco  dopth  and  rocoivor  dopth. ’ 
road(«,o)  sourco.z,  rocoivor .z 

orito(o.400) 

format  (  'Entor  transmit  boam  pattom  typo:’, 

*  1  >  Omni  boam  pattom 
’  2  ■  tablo  of  boam  pattom') 
road(«,«)  tr.boan 

Hardsiro  for  noo  tho  valuo  of  tr.boam  oqual  to  1. 
tr.boam  ■  1 

if  (tr.boam  .no.  1)  than 

orito(o,o)  'Entor  tho  namo  of  tho  filo  ohich  holds 
tho  transmit  boam  pattom. ' 
road(«.o)  filo.tr 

«rito(«,o)  'Entor  tho  anglo  for  tho  orientation  of  tho  transmittor 
(.*  is  CU  from  a  lino  from  tho  sourco  to  tho  rocoivor) .  ’ 
road(o,o)  eriontation_tr 

oriontation.tr  oriontation.tr/rad 

writo(o,o)  'Entor  tho  olovation  of  tho  transmitter  C't’  is  up).’ 
road(o,«)  olov.fix.tr 

olov.fis.tr  ■  olov.fis.tr/rad 

ondif 

vrito(«,S00) 

foraat  (  'Enter  rocoivo  boam  pattom  typo:  ’, 

’  1  ■  omni  boam  pattom 
’  2  ■  tablo  of  boam  pattern’) 
road(o,o)  roc.boam 

if  (roc.boam  .no.  1)  then 

writ«(o,o)  'Enter  tho  namo  of  tho  filo  shich  holds 
tho  rocoivo  boam  pattern. ’ 
road(o,o)  filo.roc 

*rito(o,o)  'Enter  tho  anglo  for  tho  orientation  of  tho  rocoivor 
(♦  is  CV  from  a  lino  from  tho  rocoivor  to  tho  source).’ 
road(o,o)  oriontation.roc 

oriontation.roc  ■  oriontation.roc/rad 

*rito(o,o)  'Enter  tho  olovation  of  tho  rocoivor  i*  is  up).' 
road(o,o)  olov.fix.roc 

olov.fix.roc  ■  olov.fix.roc/rad 
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277  andil 

278 

279  iritaC*,*)  ’Ent«r  th«  naa*  of  tho  fil*  ihich  holds  tho  onvolops. ’ 

280  ro4d(«,«)  fils.snv 

281  C 

282  orit«(«,*)  ’Entor  tho  noao  of  tho  filo  ohich  will  hold  tho  output.’ 

283  rood(«,o)  filo.out 

284  C 

285  «rito(*.o)  'Entor  tho  noM  of  tho  filo  which  sill  hold  tho 

286  1  scsttoring  strongth. ’ 

287  rood(o.«)  filo.ss 

288  C 

289  opon  (5,  filo  ■  filo.onv) 

290  opon  (6.  filo  ■  filo.out) 

291  opon  (7,  filo  >  filo.ss) 

292  C 

293  C  Roods  in  boon  pottoms  (if  nocossory): 

294  C 

295  if  (tr.booa  .no.  1)  thon 

296  opon  (4.  filo  ■  filo.tr) 

297  do  2  i>0.6SSl 

298  rood(4.*)  boaa_tr(S*i«^l) ,  booa_tr(S*i«’2) ,  boon  tr(S*i«3), 

299  1  booa.tr(S*i>4) .  boon.tr(Soi«S) 

300  2  continuo 

301  C 

302  rood(4,«)  booa.tr(32761) 

303  C 

304  C  if  flog  tr.boon  ••  l  thon  givo  on  onnidiroctionol  boon  pottom 

305  C 

306  olso 

307  do  S  i«l. 32761 

308  boon.trd)  ■  0.0 

309  S  continuo 

310  C 

311  oriontotion.tT  »  0.0 

312  olos.fix.tr  ■  0.0 

313  C 

314  ondif 

315  C 

316  if  (roc.boon  .no.  l)  thon 

317  opon  (3,  filo  ■  filo.roc) 

318  do  8  i-0.6SSl 

319  rood(3,o)  boon.rocCSoiol) ,  boon_roc(5*io2) ,  boan.rocCSoioS) , 

320  1  boon.roc(Soi>4),  boan.rocfSoioS) 

321  8  continuo 

322  C 

323  roodO.o)  boon.roc (32761) 

324  C 

325  C  if  flog  roc.boon  ■  1  than  giro  an  onnidiroctionol  boon  pottom 

326  C 

327  olso 

328  do  7  i-1. 32761 

329  ^on.roc(i)  •  0.0 

330  7  continuo 

331  C 

332  oriontotion.roc  >  0 

333  olos.fix.roc  >  0 

334  C 

335  ondif 

336  C 

337  C  Roods  in  tho  onsolopo  losol  fron  tho  dots: 

338  C 

339  do  10  i-0,((longth/S)  -  1) 

340  rood(5,o)  roclosoKSoiol),  roclosol(5oi42) ,  roclosol(5oi43) , 

341  1  roclosoKSoiM) ,  roclosoKSoioS) 

342  10  continuo 

343  C 

344 

345 

346  C 


orito(6.o)  sourco.x,  sourco.y 
writs (6,0)  rocoiswr.x,  rocwisor.y 
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347 

348 

349  C 

350  C 

351  C 
353 

353  C 

354 

355 

356 

357 

358 

359  C 

360  30 

361  C 
363  C 

363  C 

364 

365 

366 

367 

368 

369 

370 

371 

372  C 

373 

374  C 

375  C 

376  C 

377 

378 

379 

380 

381 
383  C 

383 

384  C 

385  C 

386  C 

387 

388 

389 

390 

391 
393 

393 

394 

395 

396 

397 

398 

399 

400  liO 

401  100 
403  C 

403  C 

404  C 

405  C 

406 

407  C 

408  C 

409  C 

410  C 

411  C 

412 

413 

414 

415 


writ«(7,«)  ’Tia«  (s«c)  Bist.  (dag)  Inc.  (dag)  Scat,  (dag) 

1  Biat.  SS  (dB)’ 

Finds  tha  allipsa  on  tha  surfaca  for  a  particular  tisa  (tisa). 

do  20  i>  start,  and.  intarval 

countar  ■  0 
araa.fina  ■  0.0 
araa.finar  -  0.0 
■ash  ■  aash.original 
flag  >  0 

if  (flag  .aq.  0)  than 

Sat  saarch  paranatars 

tina  *  i/f 

round.trip.outar  ■  iac/f 

round.trip.innar  ■  (i  -  (duration* 1000) )«c/f 
saarch.los  ■  round. trip_outar*-0. 51  -  saparation/2.0 
saarch.hi^  ■  round.trip_outar«0.51  *  saparation/2.0 
saarch.y.Ios  ■  round. trip_outar*-0. SI 
saarch.y.high  *  round.trip.outaraO.Sl 
■ash.y  ■  sash 

alsa 

Sat  saarch  paraaatars  to  around  tha  high  pt. 

saarch.los  ■  allipsa.x(location)  -  (saparatioii/6.0) 
saarch.high  *  allipsa.xdocation)  *  (saparation/6.0) 
aaarch_y.los  <■  allipsa.y(location)  -  (saparation/6.0) 
saarch.y.high  ■  allipsa.ydocation)  *  (saparation/6.0) 

■aah.y  ■  mash 

andif 

Saarch  along  x  for  y  saluas  batsaan  tha  too  allipsas 

do  100  X*  saarch.los,  saarch.bigb.  Bssh 

do  110  y>  saarch.y.los.  saarch.y.high .  Bash.y 

dist.2.snrf  •  dsqrt((x  -  sottrca.x)**2  ♦  (y  -  sourca.y)**2  ♦ 

1  (sourca.x>**2) 

dist.2.rac  ■  dsqrt((x  -  rocaisar.x)**2  *(y  -  racaiTar.y)**2  ♦ 

1  (racai*ar.x)**2) 

dist  ■  dist.2.surf  *  dist.2.rac 

if  (dist . It . round. trip.outar . and. dist . gt . round.trip.innar ) than 
countar  ■  countar  *  1 
allipsa.x(coontar)  ■  x 
allipsa.y(countar)  ■  y 
<ndif 

if  (countar  .gt.  99999)  goto  999 
continua 
continua 

Calculatas  tha  baaring  and  alasation  batsaan  tha  sourca 
and  tha  allipsa  pts. 

do  130  j*l, countar 

At  this  point,  naka  a  changa  in  coordinatas  so  that  tha  MU  of  tha 
transaittar  is  pointing  along  tha  pos.  z-axis  and  tha  transnittar 
is  locatad  at  tha  origin  (0,0,0). 

allipsa.x.tr(j)  •  ((allipsa.x(j)  ♦  saparation/2.0)* 

1  cos(oriontation.tr) 

1  -  (allipsa.y(j)*sia(oriantatioin.tr))) 

1  *cos (alav.fiz.tr)  *  (■ourca_z)*sin(alav.fix_tr) 
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416  C 

417 

418 

419  C 

420 

421 

422 

423 

424  C 

425  C 

426  C 

427 

428 

429  C 

430 

431 

432  C 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442  C 

443 

444 

445 

446  C 

447 

448 

449 

450  C 

451 

452  C 

453 

454 

455 

456  C 

457  C 

458  C 

459  C 

460  C 

461 

462 

463 

464 

465  C 

466 

467 

468  C 

469 

470 

471 
473 

473  C 

474  C 

475  C 

476 

477 

478  C 

479 

480 

481  C 

482 

483 


•llipM.y.trCj)  ■  j)  ♦  ••p4r«tioa/2.0)«sin(ori*nt»tion_tr) 

1  *  (allipa^.yC j)«cos(ori«ntation_tr)) 

•llips«_z_tr(j)  >  ((«llipa«_K(j)  *  separat ion/3.0)* 

1  coa(ori*ntation_tr) 

1  -  (*llips*_y(j)*ain(ori*ntation_tr))) 

1  •ain(*l*v_f ix.tr)  -  (aourc*_z)«coa(*l*v.lix_tr) 


Kak*  changaa  to  tho  dolta’s  to  ralfact  th*  coordinata  chango 


80urca_x.tr  *0.0 
80urca_y.tr  *0.0 

racaivar.x.tr  *  aaparation 
racaivar.y.tr  ■  0.0 


dalta.aCj) 
dalta.a( j) 
dalta.aCj) 
dalta.bCj) 
dalta.bCj) 
dalta_b(j) 
dalta.cCj) 
dalta.cCj) 
dalta.cCj) 


Caourca.x.tr  -  racaivar.x.tr) 

dalta.aCj)  ♦  Caourca.y.tr  -  racaivar_y.tr)**2 

daqrtCdalta.aC j ) ) 

C8ourca_x_tr  -  allipaa.x.trC j))**2 
dalta.bCj)  ♦  Caourca.y.tr  -  allip8a_y_trCj))**2 
daqrt  Cdalta.bC j  > ) 

Cracaivar.x.tr  -  allip88_x_tr(j))**3 
dalta.cCj)  ♦  Cracaivar.y.tr  -  allip8a_y_trCj))**2 
daqrt Cdalta.c Cj ) ) 


bmg_trCj)  ■  Cdalta_cCj))**2  -  CCdalta_aCj))**2  ♦  Cdalta_bCj))**2) 
bmg_trCj)  ■  bmg_trCj)/C-2.0*dalta.aCj)*dalta.bCj)) 
bmg.trCj)  ■  acoaCbzng_trCj))*rad 


if Callipaa.y.trCj)  'gt.  0.0)  than 
bmg.trCj)  ■  -l*bmg.trCj) 
andif 


bnig_trCj)  •  ifi»Cbmg.trCj)) 

alav.trCj)  ■>  datan2C8ourca_z,dalta.bCj))*rad 
alav.trCj)  ■  ifixCalav.trCj)) 
alav.trCj)  ■  -l**lav.trCj) 

kt  thia  point,  naka  a  changa  in  coordinataa  ao  that  tha  KU  of  tha 
racaivad  array  ia  pointing  along  tha  nag.  x-axia  and  tha  array  ia 
locatad  at  tha  origin  CO, 0,0) 

allipaa.x.racC j)  *  CCallipaa.xC j)  -  8aparation/2.0)* 

1  coaCoriantation.rac) 

1  -  Callipaa_yCj)*8inCoriantation.rac))) 

1  acoaCalav.fix.rac)  -  Cracaivar_z)*8inCalav_fix_roc) 

allipaa.y.racCj)  ■  Callipaa.xCj)  -  aaparation/2.0)*8inCoriantation_rac) 
1  *  Callipaa_yCj)*co8Coriantation_rac)) 

allipaa.z.racCj)  ■  -l*C(allipaa.xCj)  ”  aaparation/2.0)* 

1  eo8(oriantation_rac) 

1  -  Callipsa_yCj)*sia(oriaatation_ra.c))) 

1  aainCalav.fix.rac)  -  (racoivar_z)*coaCalav_fiz_rac) 

Maka  changaa  to  tha  dalta’a  to  raflact  tha  coordinata  changa 

aourca.z.rac  ■  -l*aaparation 
aourca.y.rac  ■  0.0 

racaivar.x.rac  *0.0 
racaivar.y.rac  *0.0 


dalta_a(j)  ■  Caourca.x.rac  -  racaivar.x.rac) **2 
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485 

486 

487 

488 

489 

490 

491 

492  C 

493 

494 

495 

496  C 

497 

498 

499 

500  C 

501 

502  C 

503 

504 

505 

506  C 

507  130 

508  C 

509  C 

510  C 

511 

512  C 

513  C 

514  C 

515  C 

516 

517  C 

518 

519 

520 

521 

522 

523  C 

524 

525 

526  C 

527 

528 

529  C 

530 

531 

532 

533 

534 

535 

536  C 

537 

538 

539 

540 

541 

542 

543 

544 

545  C 

546  135 

547  C 

548  C 

549  C 

550  C 

551  C 


dalta.aCj) 
d«lta_b(j) 
d«lta_b(j) 
d«lt«_b( j) 
d«lta_c( j) 
dalta.cCj) 
dalta_c(j) 


<l*lta.a(j)  ♦  (aourca.y_rac  -  racaivar.y  rac)*»2 
daqrt(dalta.a(j)) 

(aourca.a.rac  -  allipaa.s_rac(j))**2 
dalta.b(j)  ♦  (aourca.y.rac  -  allipaa.y_r«c(j))**2 
daqrt(dalta.b(j) ) 

(r*cai»ar_x,rac  -  allipa«_x_rac(J ) )«*2 
<**lta.c(j)  ♦  (racaiwar.y.rac  -  allipsa_y_r«c(j))«*2 
dsqrt (dalta.cCj) ) 


brng_rac(j)  •  (dalta.b(j))*a2  -  ((dalta_a(j))**2  ♦  (dalta.c( j))**2) 
brng.rac(j)  -  bmg.rac(j>/(-2.0*dalta.a(j)*dalta_c(j)) 
brng.rac(j)  •  acoa(brng_rac(j))*rad 


if (allipsa.y_rac(j)  .It.  0.0)  than 
brag_rac(j)  ■  -l*bmg_r«c( j) 
andil 


brng.rac(j)  •  ifix(brng_rac(j)) 

alav.rac(j)  •  datan2(allipaa.z_rac(j) ,dalta_c( j) )*rad 
alav.rac(j)  ■  ifix(alav_rac(j)) 
alav.racCj)  ■  -l*alaw_rac( j) 

continua 

Raaat  aax  aach  tiaa  through  tha  loop. 

■ax  >  >1000.0 

Saarch  through  tha  points  on  tha  allipsa  tor  tha  loaast  loss  (  tha 
sua  of  tha  transaission  lossas  and  tha  baaa  gains). 

do  135  jal, countar 

if (bmg_tr<j)  .gt.90  .or.  bmg.tr(j) .lt.-90 
.or.  bmg.rac(j). gt.90  .or.  bmg.rac(j) .lt.-90)  than 
tast(j)  -  -1000.0 
go  to  135 
andif 

pointar  ■  (bmg_tr(j)a90)ai81 
pointar  -  pointar  *  alav.tr(j)  91 

pointar.rac  ■  (bmg_rac(j)^90>»18l 
pointar.rac  ■  pointar.rac  ♦  alav.racCj)  +91 

tastCj)  ■  baaa.tr (pointar) 

tastCj)  ■  tastCj)  ♦  baaa.rac (pointar.rac) 

tastCj)  ■  tast(j)  -  iOadloglO(dalta.b(j)aa2  ♦  (■ourca.z)*a2) 

tastCj)  ■  tastCj)  -  10adlogl0(dalta.c(j)a*2  ♦  (racaivar.z)*a2) 

tastCj)  ■  tastCj)  -  sttn*dsqrt(dalta.b(j)**2  ♦  (sourca.z)a*2) 

tastCj)  ■  tastCj)  -  attnadsqrtCdalta.cC j)*a2  ♦  (racaivar.z)««2) 

if  (tastCj) .gt.aax)  than 
location  ■  j 
■ax  -  tast(j) 
baaring.hit  ■  bmg.tr(j) 
alavation.hit  ■  alav.trCj) 
position.tr  ■  pointar 
position.rac  ■  pointar.rac 
andif 

continua 

Chack  to  saa  if  flag  is  zaro,  if  so,  sat  aqual  to  1.  calculata 
tha  arts  and  drop  tha  aash  sizs.  If  flag  is  nonzaro,  calculata 
tha  araa  and  coapara  with  coarsar  aash  slza  araa. 


553 

553 

554  C 

555  C 

556  C 

557  C 
556 

559 

560 

561 

563  180 

563  C 

564 

565 

566 

567  C 

568 

569  C 

570  C 

571  C 
573  C 

573 

574 

575 

576 

577  300 

578  C 

579  C 

580  C 

581  C 
583  C 

583  C 

584  C 

585 

586 

587 

588 

589 

590 

591 
593 

593 

594 

595  C 

596 

597  C 

598  C 

599  C 

600  330 

601  C 
603  C 

603  C 

604 

605 

606 

607 

608 

609 

610 
611 
613 

613 

614 

615  C 

616 

617 

618 
619 
630 


if  (flag  .aq.  0)  than 
flag  ■  1 

Count  all  pointa  on  tha  aurfaca  that  ara  aithin  3  dB  of 
■ax  and  calculata  tha  araa. 

do  180  j  ■  1 .  countar 

if  ((nax  -  3.0)  .It.  taat(j))  than 
araa.fina  -  araa.fina  *  (■aah)*a3 
andif 
continua 

■ash  ■  ■aah/3.0 

countar  •  0 
goto  30 

alaa 

Count  all  points  on  tha  aurfaca  that  ara  aithin  3  dB  of 
■ax  and  calculata  tha  araa. 

do  300  j  ■  1.  countar 

if  ((nax  -  3.0)  .It.  tast(j))  than 
araa.finar  ■  araa.finar  *  (■aBh)a*3 
andif 
cont inua 

Aftar  tha  araa  has  baan  calculatad,  conpara  aith  tha  araa  calculatad 
aith  tha  prior  sash  sixa  (pass  if  aithin  */-  lOX) ■  This  naans  that 
tha  first  accaptabla  araa  can  ba  at  nash.original/S . 0 .  If  tha  araas 
ara  not  aithin  tha  giasn  tolaranca.  drop  tha  nash  siza  by  a  factor  of 
tao  and  do  it  again. 

if  ( araa.f inara l . 3S . gt . araa.f ina . and . araa.f inaraO . 8 . 1 t . araa.fina ) than 
araa  *  araa.finar 

goto  330 

also 

araa.fina  ••  araa.finar 
araa.finar  ■  0.0 

■ash  ■  ■ash/3.0 

countar  ■  0 
goto  30 
andif 

andif 

Calculata  tha  bistatic  surfaca  scattaring  strangth  and  srita  it. 
scat.str  ■  raclaToKi-bagin)  -tast (location)  -trlaval  -lOadloglO(araa) 
Calculata  all  tha  anglos  and  «rita  than  out. 
dalta.adocation)  ■  (sourca.x  -  racaivar.x)  **3 

dalta.a(location)  «  dalta.adocation)  a  (sourca.y  -  racaivar.y)aa3 
dalta.adocation)  *  dsqrt(dalta.a(location)> 
dalta.bdocation)  ■  (sourca.x  -  allipsa.xdocation) )«*3 
dalta.bdocation)  ■  dalta.bdocation)  * 

1  (sourca.y  -  allipaa.y(location))aa3 

dalta.bdocation)  •  dsqrt(dalta.b(location)) 

dalta.c (location)  ■  (racaivar.x  -  allipsa.xdocation) )*a3 

dalta.c (location)  ■  dalta.c (location)  * 

1  (racaivar.y  -  allipsa.y(location))**3 

dalta.cdocation)  -  dsqrt (dalta.c (location)) 

bist.angla  ■  (dalta.adocation)  )*«3 
bist. angle  ■  bist.angla  -  ((dalta.bdocation)) ••3 
1  *  (dalta.cdocation)  )«a2) 

bist.angla  ■  bist.angla/(~3.0*dalta.b(location)adalta.c(location)) 
bist.angla  -  acos (bist.angla )arad 
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631 
623  C 
623 

634  C 

635 
626  C 
62T 

638  260 

639  C 

630 

631  C 
633 

633  C 

634  999 

635  C 

636  20 

637  C 

638 

639 


bist.angl*  ■  180.0  -  dab«(bist_angl«) 

th«t>.n4ught  ■  clatui2(sourc«_z.d«lta.b(location))*r*d 

th«ta  -  datuk2(recciv«r.z.d«lta_c(location))*rad 

■rit«(7,260)  tiaa,  bist.angla,  thata.naught.  acat.str 

foraat(16.3.4z,lf6.i.6z.lf6.1,7z.lf6.1.7z,f6.1) 

arit«(6.«)  •llipa«_z(location) .  allipaa.y(location) ,  tiaa 

goto  20 

arita(7,«)  'At  atop.  Countar  *  countor 

cont inua 

atop 

and 
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